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November 1988 ^ough September 1992. The project was funded by the National 
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Program (In-STEP), and the contract was administered by NASA Lewis Research 
Center. 


NASA P*roject Manager: 

Richard H. Knoll 

Key NASA Personnel: 

Todd Torn 
John Aydelott 
Dr. M. Mojibul Hasan 
Dr. Chin-Shun Lin 


Boeing Program Manager: 

Dr. Jere S. Meserole 

Boeing Principal Investigator: 

Michael D. Bentz 

Key Boeing Personnel: 

Dan Luteran 
Glenn Sikkema 
Donald Skene 
Gary Frankfurt 
Gary Brown 


Deputy Project Manager 
Science Advisor 
Project Scientist 
Project Scientist 


mechanical design 
electronics design 
technician 
technician 
software 


Data are presented in System Internationale (SI) units, with English units also given 
where appropriate. 

The complete set of video imagery obtained on this effort is available as an adjunct to 
this final repOTt through the NASA Lewis Research Center videotape hbrary. A request 
form for this video supplement is shown on the last page. 




CONTENTS 


FOREWORD i 

CONTENTS iii 

SYMBOLS V 

SUMMARY 1 

INTRODUCTION 3 

BACKGROUND 7 

OBJECTIVES 13 

EXPERIMENT APPROACH AND REQUIREMENTS 15 

General 15 

Instrumentation 18 

EXPERIMENT APPARATUS AND OPERATION 21 

STS Accommodations 21 

Structure 24 

Fluid System 27 

Power Supply 32 

Controls and Data Storage 33 

Instrumentation 34 

Photography 35 

Operating Timeline 37 

LOW-G AIRCRAFT TEST 39 

Aircraft Test Objectives and Approach 39 

Aircraft Test Results and Conclusions 40 

IN-SPACE TEST 43 

Payload Performance 43 

Data Summary 47 

ANALYSIS AND DISCUSSION 63 

Flow Patterns 63 

Mixing Times 66 

Comparison With Mixing Models 74 

Heater Plate Convection 80 

Heating Hiase Effects 81 

CONCLUSIONS 83 

REFERENCES 85 

APPENDIX A A-1 


m 


( \ 




PRECEDING PAGE BLANK NOT FILMED 




SYMBOLS 


A liquidAapor interface area 

Boj Bond number of jet at liquid/vapor interface, pgDj2/a 

Cpi liquid specific heat 

Dj jet diameter at liquid/vapor interface 

Do jet diameter at nozzle exit (nozzle inside diameter) 

Dtank tank diameter 

F flow characterization parameter (see equation 6) 

g acceleration 

gpm gallons per minute 
h convective heat transfer coefficient 

hfe jet submergence (distance fiom nozzle exit to liquid/vapor interface) 

ID inside diameter 

ki liquid thermal conductivity 

1 liters 

lAnin liters per minute 

m mass 

N dimensionless buoyancy parameter 

Nu Nusselt number based on heater length 
P tank puessure 

Pnc sum of partial pressures of noncondensible gases 

Pri liquid Prandtl number, CpiPi/ki 

Q volumetric flow rate 

Rej jet Reynolds number, VqDo/Vi 

T fluid temperature 

Tsat saturation temperature 

AT maximum fluid temperature spread 

t time 

tmix time required to mix fluid (see pages 5-6 for specific definitions of mixing times) 
V ullage volume 

Viiquid volume of bulk liquid 
Vo jet velocity at nozzle exit 
W watts 

Wej Weber number of jet at liquid/vapor interface, piVo^Ro^/oDj 

Zb liquid level (distance from liquid/vapor interface to mixer-end of tank) 

“k latent heat of evaporation 

Pi liquid dynamic viscosity 

Vi liquid kinematic viscosity 

p ullage density 

Pi liquid density 

o siuface tension 

X time constant for exponential decay of temperature and pressure 
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SUMMARY 


The Tank Pressure Control Experiment (TPCE) is a space experiment designed to 
develop an essential cryogenic fluid management technology: control of storage tank 
pressures in the absence of gravity by forc^-convective mixing. Axial-jet recirculation 
mixers are candidates for controlUng pressure rise in cryogenic fluids stored for future 
long-duration noissions, such as in an orbiting propellant depot. Positive control of tank 
pressure under all fluid orientations and conditions is critical, but the addition of 
excessive energy to the fluid from such devices must be minimized to avoid vehicle 
disturbances and ciyogen boiloff losses. 

The objectives of TPCE were to characterize the fluid dynamics of axial-jet-induced 
mixing in low gravity, to evaluate the validity of empirical mixing models, and to 
provide data for use in developing and validating computational fluid dynamic models of 
mixing processes. The experiment used refrigerant- 113 as the test fluid, at near- 
saturation conditions, with the equilibrium pressure approximately 4.2 kPa (0.6 psi) 
above the vapor pressure. This simulates the fluid dynamic and thermodynamic 
behavior of subciitical cryogenic fluids in space that are stored at or near their boiling 
points. The 13.7-liter, 25-cm-diameter tank was filled to an 83% fill fiaction. 

Immersed heaters created pressure increases and thermal gradients. A 1.0-cm-diameter 
nozzle produced an axial jet that was mmed on to mix the fluid and return the pressxue to 
a near-equilibrium level. Video cameras recorded the liquid/vapor interface shapes and 
flow patterns produced as a function of mixer flow rate, while temperatures and 
pressures were recorded to assess the mixing and pressure control effectiveness. The 
experiment was packaged as a Complex Autonomous Payload utilizing the Get-Away 
Special container in the Shuttle payload bay as its carrier. 

TPCE accomplished aU of its objectives in flight on the Space Transportation System 
(STS) in August 1991. Two of the flow patterns observed generally agreed with prior 
drop tower test observations and correlations. The jet penetrated the ullage and formed a 
symmetrical recirculating flow pattern at jet Weber numbers of ^proximately 4.8 and 
above. Weber numbers of 1.4 or less pr^uced a symmetrical, nonpenetrating pattern. 
Mixing was found to be effective at both of these flow patterns, in terms of the ability to 
produce pressure collapse and thermal equilibration. The penetrating jet provided the 
most rapid pressure reduction and equilibration. Jet Weber numbers between 1 .4 and 
4.8 resulted in asymmetric flows which were less consistent in their mixing 
effectiveness. A closed-form equation derived from a simple thermodynamic model was 
found to provide an approximate prediction of the pressure reduction time as a function 
of mixer parameters, tank size, and fluid thermophysical properties. 

A low-energy axial jet mixer was found to be effective and consistent at mixing, 
enhancing heat and mass transfer between the phases, and controlling tank p>ressure. 
Active mixing — which offers increased predictability and controllability in space 
cryogenic systems compared with passive pressure control designs — was demonstrated 
to be possible using a low kinetic energy jet that incurs no significant boiloff penalty. 




INTRODUCTION 


Reliable control of tank pressure is essential to the storage and handling of cryogenic 
fluids in space. A number of candidate future missions will require the use of large 
quantities of cryogenic fluids and advanced concepts for storing and handling them. 
Applications such as reusable space transfer vehicles, manned Mars stages, and orbital 
propeUant depots will require the storage of cryogenic fluids for greater lengths of time 
than has been demonstrated to date. 

Any tank of cryogenic liquid stored in space is subject to heat addition from its 
environment, often at varying rates and with a nonuniform distribution around the tank. 
These heat inputs raise the tank pressure by evaporating liquid and/or warming the 
vapor. The absence of strong gravity-driven natural convection can lead to nonuniform 
distribution of the heat in the liquid. As the pressure in a closed tank rises in response to 
the heating, its rate of increase depends strongly on the distribution of the heat within the 
fluid. A nonuniform temperature distribution results in higher pressure than would be 
found with an isothermal fluid. Localized "hot spots" in the fluid drive the pressure up 
faster and to higher levels because the tank pressure is strongly driven by the vapor 
pressure of the warmest liquid in the tank. If the heat input is instead distributed 
uniformly throughout the fluid, then the pressure rise is ininimized. The pressure rise 
rate in a tank of well-mixed cryogen can be more than an order of magnitude lower than 
that in a tank of ciyogen that is not in equilibrium. 

The near-absence of gravity not only reduces natural convection, but also makes it 
difficult to relieve the resulting pressure increase by the conventional venting of vapor 
since the locations of the liquid and vapor are not generally known with certainty. For 
long-term storage, pressure must therefore be controlled by thermodynamic means such 
as a refrigerator or a thermodynamic vent system (TVS). A TVS, which is essentially 
an open-loop refrigeration cycle, j^ovides cooling by expanding a small amount of 
liquid or vapor through an expansion valve, passing it through a heat exchanger to 
absorb heat from the remaining fluid, and then venting the warmed vapor overboard. 

The effectiveness of either a TVS or a refrigerator depends on the distribution of cooling 
to all regions of the tank. To accomplish this, cither system can be coupled with a mixer 
using a compact heat exchanger that is integral to the mixer device, as illustrated 
schematically in figure 1. An internal heat exchanger employing forced craivection 
offers lower mass than a passive system using only an external TVS-cooled shield that 
relies on conduction and natural convection. 

For short-term storage periods (such as the several days that might elapse between the 
launching of an upper stage to orbit and the firing of its engines), mixing alone could be 
an adequate means of controlling pressures. A relatively small amount of induced liquid 
motion may be sufficient to break up the hot spots and provide forced-convection heat 
transfer at a rate many times higher than that for free convection and conduction in low 
gravity. 
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Figure 1: Pressure Control System Using Mixer 
The benefits of active mixing to space cryogenic systems are that it: 

1 . Reduces pressure-rise-rate in a nonvented tank. This allows systems to be designed 
that can store cryogen for longer periods without venting or refrigeration. 

2 . Enables the use of compact heat exchangers. Forced convective mixing increases 
the effectiveness of a refrigerator or TVS and allows the use of an intanal compact 
heat exchanger, thus reducing weight coir 5 )arcd with a passive system. 

3 . Can serve as a backup for passive cooling systems. If cooled extranal shields are 
used by a TVS or refrigerator to intercept heat leaks, rather than an internal heat 
exchanger, the ability of such a system to cool the fluid depends primarily on heat 
conduction. An active mixer can be added as a backup to this passive system to 
increase confidence in its ability to provide positive control of tank pressure. 

4. May reduce on-orbit fluid tranter time. The rate of a non-vented fill depends on the 
rate at which vapor can be condensed, which in turn depends on the rate of 
dissipation of the heat of compression and condensation away from the liquid/vapor 
interface. Although the incoming liquid stream will induce some circulation, under 
some circumstances additional forced circulation could speed the condensation and 
thereby reduce the fill time or increase the final liquid fill level that is attainable. 

5 . Ensures delivery of uniform-temperature liquid to engines. For propellant tanks, 
mixing can provide a more uniform propellant supply to the engines, which 
potentially reduces pressurization requirements and residual propellant mass. 

6. Makes the fluid state predictable. Periodic mixing reduces or eliminates the potential 
for sudden pressure rise or drop due to nonequilibrium fluid states when the fluid is 
subjected to sudden motions. 

The problem is not whether mixing systems can be built that provide adequate 
circulation, but how to mix efficiently as well as reliably. Mixers have been designed 
for use in previous cryogenic space systems (the Saturn S-IVB stage and the cancelled 
Shuttle/Centaur stage [ref. 1]) based on extensive ground testing. These, however, had 
to be designed conservatively to compensate for the lack of certainty about low-g mixing 
behavior. 
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One penalty for conservatism is that a high-velocity mixing jet could result in fluid 
frothing and in undesirable vehicle disturbances, particularly for applications in which 
propellant is a large percentage of the vehicle mass. Another penalty is that the ener^ 
used by the device eventually dissipates as heat in the fluid, resulting in an increase in 
propellant boiloff. Thus there is a need to minimize the kinetic energy added by the 
mixer. 

For example, the Shuttle/Centaur hydrogen tank mixer is rated at 15W total power input 
while operating (ref. 1). Of this 15W, the jet kinetic energy rate is approximately 0.4W 
and the rest is due to pressme losses in the TVS, mixer inefficiency, and motor losses. 
In long-term storage applications, a mixer dissipating even 0.4W would represent a 
significant heat addition, resulting in 25 kg (55 Ibm) per year additional boiloff of 
hydrogen or 60 kg (132 Ibm) per year of oxygen. Previous mixer performance models 
indicate that this energy input could be reduc^ at least one order of magnitude by using 
a larger jet diameter and lower jet velocity. Lx)w-kinetic-energy jet mixing, however, 
was unproven in its effectiveness with a low-g flitid orientation prior to the present 
effort. Without adequate low-g data, mixers would continue to be designed 
conservatively (i.e., with higher-velocity flow) resulting in both mass and power 
penalties. 

Many investigators over the past twenty-five years have recognized the need for 
understanding mixing in space-vehicle tanks. Their studies have included small- and 
large-scale thermal tests in one g with a number of test fluids, plus low-g mixing 
simulations and drop tower tests. Several analytical models and empirical correlations 
were developed, as described in the following background section, and valuable 
conclusions and understandings emerged from these efforts. However, because of the 
complexity of the fluid dynamic problem and the difficulty of studymg it in a normal- 
gravity environment, there has b^n no validated model of low-g mixing behavior. 

The following sections describe the objectives of the TPCE space experiment, the 
technical approach and requirements for meeting these objectives, the experimental 
apparatus and operating s^uence, the results of preliminary low-g testing on an aira-aft, 
and finally the results of a 25.3 hour test on board the Space Shuttle Atlantis (mission 
STS-43). The observed flow patterns and the measured mixing effectiveness are 
discussed and compared with the prior correlations, and a simplified thermodynamic 
model of pressure control effectiveness is derived and compared with the flight data. 

The complete set of digital data is presented in graphical form in Appendix A, and 
instructions for obtaining a copy of the video d^ are given. 
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BACKGROUND 


In 1968 Poth et al. (ref. 2) evaluated a variety of mixing concepts for low-g cryogenic 
pressure control devices. Concepts such as paddlewheels, dielectrophoresis, radial jets, 
axial jets, and other devices were studied. The desired capabilities of a mixer includ^: 

1 . The ability to provide thorough mixing with minimal kinetic energy addition. 

2 . The ability to disrupt (penetrate) the ullage, to reach all regions of the tank 
with the mixer flow, and to sweep the tsmk walls. 

3 . The ability to provide effective mixing at all fill levels. 

4. The avoidance of swirling motions and non-symmetrical forces exerted on 
the tank walls, to avoid disturbances to the vehicle. 

5 . The ability to operate intermittently (non-continuous operation) if desired. 

6. Compatibility with liquid hydrogen. 

The mixing device also should be compact, lightweight, and capable of being integrated 
with a compact heat exchanger. Based on all of these evaluation criteria, the axial liquid 
jet was found by Poth to be most suitable for spacecraft propellant tank pressure control. 

Most prior work to predict low-g mixer performance has been based on analysis and on 
normal-gravity experiments. An industrial mixing equation, published in 1948 by 
Fossett and Prosser (ref. 3), uses a nondimensional correlation of the time required to 
blend liquids: 

tmix^pDo _ g (j) 


Here tmix is defined as the time required to distribute a chemical species throughout a 
liquid. Vo and Do are the velocity and diameter of the jet, respectively, at the outlet of 
the mixer nozzle, and Dtank is the diameter of the vessel. A later correlation by Fox and 
Gex (ref. 4) incorporates more parameters including the kinematic viscosity of the 
liquid, V, and the liquid level, Zb. Their correlation can be rearranged in terms of the 
same dimensionless grouping for comparison: 

tmixVoDo ^ (Re>2,000) (2) 

^tank 

The dependence on viscosity and velocity in the right-hand side of this equation is weak, 
and if Zb, Do, and Dtank scale together then the dependence on tank size is also weak. 
This equation, however, can lead to predicted mixing times for large liquid hydrogen 
tanks that are an order of magnitude greater than those predicted by equation (1). 

Whereas equations (1) and (2) were derived for the mixing of various industrial 
chemical compounds, Poth and Van Hook (refs. 5, 6) conducted a series of axial jet 
stratification-reduction tests with heated water in tanks 0.3 m (1 ft) and 3 m (10 ft) in 
diameter. The same dimensionless grouping was used to present mixing time, defined 
as the time required for the axial jet mixer to bring all liquid temperatures measured to 
within ten percent of each other (results were also given for 20% and 5% convergence 
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times, as well as pressure collapse times). The results of these normal gravity tests were 
presented in terms of a dimensionless buoyancy parameter N: 

'■^7°°° = /(N) (3) 


As the buoyancy parameter N (which accounts for the thickness of the stratified layer as 
well as its buoyancy) approaches zero, the functiony(N) approaches a constant, having 
a value between approximately 1 and 4 (for 10% temperature convergence). 

WoUen (ref. 7) tested jet mixing in simulated low gravity by obsCTving the time required 
for an acid/base mixture to become homogeneous enou^ for a pH-indicator dye to 
become colorless. He reported a constant value of the same dimensionless mixing time, 
for sufficiently high Reynolds numbers: 


hnixVoDo 



3.46 


(Re > 5,000) 


(4) 


Meserole et. al (ref. 8) measured thermal mixing times in refiigerant-1 1 in normal 
gravity in a 2S-cm diameter tank with varying nozzle diameters. At flow rates high 
enough to overcome the retarding effect of buoyancy (i.e. at Reynolds numbers greater 
than 15,000) the dimensionless mixing time agreed with the value reported by Wollen. 


Dominick (ref. 9) tested a side-mounted jet in refrigerants 11 and 113 in one-g, 
measuring steady-state condensation rates at the liquid/vapor interface rather dim 
transient pressure reduction rates. He found that the submerged jet, when impinging on 
the surface, increased ullage condensation rates by a factor of seven over those obtained 
with no mixing in normal gravity. 


Bullard (ref. 10) and Sterbentz (ref. 1 1) successfully demonstrated the effectiveness of a 
compact TVS combined with an axial mixing jet in ccmtrolling the pressure of liquid 
hydrogen in normal gravity. These tests, using 1.04-m (41-in) and 2.8-m (1 10-in) 
tanks, also confirmed that Ae axial jet is more effective than a side-noounted jet. 

Pressure decay rates were controlled by the rate of heat transfer in the heat exchanger, 
and were approximately predicted by a simple one-lump thermodynamic model. 

The above experiments had a one-g liquid/vapor configuration; that is, a flat interface. 
Mixing is inherently different in low gravity because of the liquid/vapor orientation. 
Under conditions having sufficient orbital drag with minima] disturbances, such as 
large-diameter tanks at low orbital altitudes (e.g., a propellant depot), the liquid will tend 
to be settled at one region of the tank, but its interface will be highly curved by surface 
tension forces. Under other conditions the liquid position will be highly variable and 
affected by acceleration noise ("g-jitter") and boiling. Experiments using heated (non- 
isothermal) fluids in one-g are also strongly affected by buoyancy, which stratifies the 
fluid and retards mixing. 


In an early low-gravity mixing experiment Berenyi, et al. (ref. 12) conducted drop- 
tower tests of jet mixing in 10-cm (4-m) and 20-cm (8-in) spherical tanks. The objective 
of these experiments was to determine the jet velocity required to penetrate all bubbles 
and provide complete circulation in the tatik. 


Aydelott (refs. 13, 14, 15) tested an axial jet mixer in several 10-cm (4-in) diameter 
tanks of ethanol during fieefall in the NASA Lewis Research Center 5-second drop 
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tower facility. I^e was added to the jet and the mixing time, defined as the time 
required to distribute the dye throughout the bulk liquid, was determined visually. 
Mixing times were correlated by: 


= (0.09) + (0.01) ln(F) (5) 

where Q is the volumetric flow rate, Viiquid is the volume of the bulk liquid, and F is a 
flow characterization parameter defined (for completely turbulent jets) as 


^ _ (1.6)Wej - (0.5) 
^ ■ (1) + (U.6)Boj 


(Re > 1,250) 


( 6 ) 


Aydelott defines the jet Weber number, Wej, and Bond number, Boj, as the values at the 
zone where the jet meets the liquid/vapor interface, using the diameter of the jet at the 
interface, Dj, as the characteristic dimension: 


Wej = 


Pi Vq^Rq^ 
aDi 


Boj=pLue£ 

a 


(7) 

( 8 ) 


The jet Weber number defined here is the ratio of jet momentum to the surface tension 
force that resists jet penetration at the liquid/vapor interface. The jet Bond number is the 
ratio of the force of gravity or acceleration on the jet geyser to that of surface tension. 

The left-hand side of equation (5) is a different nondimensionalization of the mixing time 
than in the first four equations. This correlation was not intended to encompass a tank- 
to-jet diameter ratio other than 25: 1 and, in comparison with other test data, would not 
be expected to properly reflect the effect of jet diameter on mixing time. However for a 
constant ratio of tank-to-nozzle diameter and for a given tank shtqie and fill level the 
parameter is essentially proportional to the dimensionless mixing time parameter shown 
in equations (1) ^ough (4). This supports the use of (tmixVoDo^tai^ as a correlator 
of fluid mixing time, at sufficiently-high Reynolds numbers, for a wide range of fluids 
and tank sizes. 

It may be observed that this parameter, (tmixVoDo/D^tank). is related to the square root 
of the jet momentum rate, rather than to the kinetic energy rate, 

(’^/gjpVp^Do^. This would in^cate that mixing effectiveness is determined by the rate 
of addition of momentum (the jet's thrust), rather than by the rate of addition of kinetic 
energy (the jet’s power). Therefore it appears to be most efficient to design a jet mixer 
with minimum kinetic energy but with sufficient momentum, or (VqDo) 2, to produce the 
desired pressure control effectiveness. This would imply that the velocity of the jet 
should be minimized while the diameter, and therefore the volumetric flow rate, should 
be increased accordingly. (In designing a system, the diminishing return as the diameter 
is increased would be traded off against the size and mass of the device.) 

Aydelott (ref. 13-15) identified four jet flow patterns from his tests. Low-g Flow 
Pattern I is defined as dissipation of the jet in the bulk liquid region without breaking the 
interface (this is the regime most similar to one-g experiments). In Flow Pattern II tiie 
jet penetrates and forms a geyser, and in Pattern m tiie geyser reaches the opposite end 
of Ae tank and collects there. Flow Pattern TV is defined as flow of the jet over the 
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opposite end of the tank and along the walls back to the bulk liquid region, giving 
complete circulation throughout the tank. Pattern I was observed at values of F less than 
0.4 (the transition between patterns I and II ranges from 0.4 < F < 1.0). How Pattern 
IV was found at F > 2.5. The correlation of his mixing time paranKter with the flow 
parameter given in equation (5) was valid only fcM: Patterns I and IV. Patterns H and HI 
were observed to be irregular and less effective at mixing and therefore are considered to 
be undesirable flow conditions (although this observation is based on very short- 
duration tests). 

Since the Weber number defined in equation (7) is proportional to the jet momentum, 
normalized by dividing it by a surface tension force, Wei may be a useful parameter for 
predicting both flow pattern and mixing times. Aydelotrs flow charactaization 
parameter F successfully carrclated data at var^g g-levels including one g, but it is 
essentially proportional to Wej for cases in which the jet Bond number, B^, is much 
less than one (as will usually be the case in low g). Since the use of equations (5) and 
(6) causes difficulties at low values of Wej which give negative values of F, one can 
equivalently use Wej to correlate flow patterns at sufficiently small Ba values. 
Considering only the subset of Aydelott’s data in which Boj is small. Flow Pattern I and 
the borderline Pattern I/H are seen at Wej values between 0.31 and 1.11, Patterns II and 
in are observed at 1 . 1 6 < Wej < 2.04, and Pattern IV is seen at Wej values of 1 .33 to 
ovCT 200 (note that there is some overlap between the Pattern IE and Pattern IV ranges). 

Sterbentz and Bullard (ref. 11) give an analytic prediction of the jet inomentip required 
to penetrate the ullage, using an assumed 9* jet spread half-angle. This prediction was 
reported to agree approximately with Berenyi's (ref. 12) drop tower results. The 
penetration criterion, which is based on a simple force balance at the meniscus, can be 
rearranged in terms of Wej as: 

Wej, critical =1 

Figure 2 depicts the zone of interaction between the jet and the liquid/vapor interface. 
Poth (ref. 2) gives a more detailed jet penetration analysis that accounts for the effect of 
eddy viscosity on the jet velocity proHe. The resulting criteritMi is: 

> (1.055) hb (10) 

ct gc 

which when rearranged in terms of Wej is: 

Wej, critical ~ (0.26) (1 1) 

If one assumes, as did Aydelott, that a completely turbulent jet (Re > 1250) spreads at 
approximately 7* for the first 6.2 outlet diameters and thereafter at 1 1°, the Poth criterion 
reduces, for (Do/hb) « 1, to: 

Wej. critical = 0.69 (12) 

Thus this analytical prediction of the jet momentum required to penetrate the ullage is in 
agreement with the low-g experimental data of Aydelott (ref. 15). 

Computational fluid dynamics modelling of submerged axial jets with free surfaces has 
more recently been performed by Hochstein etal. (ref. 16), Wendl (ref. 17), and 
Wendl et al. (ref. 18). Predictions of the critical jet velocity required to penetrate the 
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ullage, made for the TPCE geometry under a subcontract to the present effort (ref. 17), 
a^ed weU with the predictions of Aydelott (ref. 15). The numerical codes, which are 
still under development, could ultimately predict thr^-dimensional flow patterns and the 
thermodynamic response to mixing (i.e., pressure and temperature) for a tank of 
arbitrary size and with any cryogenic fluid Further code development and validation 
will require low-g fluid dynamic and thermodynamic data for comparison with the 
predictions. 
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OBJECTIVES 


The experimental work described in the preceeding section was subject to the inherent 
limitations of ground-based testing. Gravity constrained the liquid/vapor geometry and 
caused buoyancy effects in the one-g tests. The low-g drop tower tests were limited to 
small size tanks and very short test durations. 

The next step toward a complete understanding of pressure control in low g is to test 
axial jet mixing of a non-isothermal fluid in a realistic low-g environment. To be 
realistic, the gravity level of the test environment must be low enough that the fluid is 
dominated by capiUaiy and inertial forces rather than acceleration forces, and the 
duration must be sufficient to observe both the fully-developed flow patterns and the 
complete pressure collapse. Except with extremely small scale test tanks, these 
conditions can only be met by an in-space experiment 

The objectives of the TPCE project are to measure the effectiveness of an axial jet mixer 
at pressure control in low g, to evaluate existing mixing models, and to provide data for 
developing and validating future models. Specifically, the experiment objective is to 
determine the effectiveness of jet mixing as a function of flow rate, where effectiveness 
is defined by; 

a. the ability of the jet to disrupt the vapor and reach all regions of the tank, 

b. the rate at which mixing equilibrates temperatures in the fluid, and 

c. the rate at which mixing reduces the pressure in the tank. 
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EXPERIMENT APPROACH AND REQUIREMENTS 


General 

The Tank Pressure Control Experiment is designed to be carried as a Shuttle secondary 
payload within a Get-Away Special (GAS) canister. This constrains the experiment to 
the use of a noncryogenic simulant fluid. Refrigerant- 1 13 (trichlorotrifluoroethane) was 
chosen because of its benign physical characteristics and because of its moderate vapor 
pressure over the range of temperatures expected during flight In addition to the higher 
boiling point, the thermophysical properties of this test fluid differ from those of 
cryogenic propellants in several important ways, particularly in the liquid's thermal 
diffusivity and Prandtl number. Also, the scale of the exp)mment and the scope of the 
project do not permit a themaodynamically-scaled exp>eriment For instance, the 
important ratio of tank thennal mass to volume cannot be held constant, nor can a 
redistic tank thermal environment be reproduced without compromising the ability to 
view the fluid through the tank walls. Therefore the emphasis of TPCE is on the fluid 
dynamics. An imderstanding of the fluid dynamic behavior is of primary impx>rtance. 
Any model of pressure control processes will be strongly driven by the predicted flow 
patterns and liquid/vapior interface areas, since these control the rates of heat and mass 
transfer. However, thermodynamic data (temperatures and pressure) are acquired in 
TPCE to allow an assessment of relative mixing effectiveness and to obtain ^ta for 
comparison with models. 

To simulate cryogenic propiellants that are stored at or near their boiling pxiint, the test 
fluid is also near its saturation pnint As a goal the total tank pressure (at equilibrium) is 
equal to the vapxir pressure of the mixed fluid, plus a minimum excess pressure which is 
required to prevent pump cavitation in the absence of gravitational head pressure. This 
excess pressure is referr^ to as the total partial pressure of noncondensible gases, or 
Pnc- Complete cavitation suppression requires a Pnc of approximately 7 kPa (1 psi) for 
TPCE's fluid system, as determined by developmental testing. To provide this amount 
with miiumiun effect on the ullage thermal capacitance it was desir^ that the Pnc consist 
primarily of helium gas. Thus the refrigerant was to be as fioe of other dissolved gases 
as possible, with approximately 7-14 kPa (1-2 psi) partial pressure of helium. (As will 
be discussed in the Payload Performance section, the actual composition of Pnc in flight 
was somewhat different than this goal.) 

The GAS container allows a total payload volume of 0.14 m^ (5 ft^) and mass of 91 kg 
(200 Ibm). TTiis constrains the geometry of the expjeriment’s fluid tank. The geometry 
of the tank was chosen based on a review of the propiellant tank dimensions proposed in 
several smdies of space transfer vehicles and proprellant depots. A tank geometry was 
chosen that is representative of liquid hydrogen tanks in these applications. This is a 
cylinder with hemispherical ends, having a length-to-diameter QJD) ratio of 1.4. 
(Typical hydrogen tank L/D’s range from 1.0 to 2.5 in the literature, while oxygen tank 
L/D's range from 0.7 to 1.0.) These considerations resulted in a 13.7-liter (0.485-ft^) 
tank with a 25.4-cm (10-in) inside diameter, as shown in figure 3. This is 
approximately a 5% to 7% scale model (by linear dimension) of a typical space transfer 
vehicle hydrogen tank. 

A jet nozzle is positioned along the tank's major axis near one end of the tank. The 
nozzle outlet is displaced 6.3 cm (2.5 in) from the end of the tank to repjresent the 
geometry of a combined mixer/TVS in an application. The mixer pumps are located 
outside of the tank due to the difficulties of developing a scale-model internal pump that 
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meets STS safety requirements. These external pumps are supplied \vith liquid from the 
tank by a liquid acquisition device (LAD) so that the jet will not contain vapor. 



Figure 3: TPCE Test Tank Geometry 

The jet nozzle is a straight tube with a 1-cm (0.40-in) inside diameter. This diameter 
was chosen so that the tank-to-jet diameter ratio is 25, for similarity to several prior 
experiments. Aydelott's tests (refs. 13 and 15) used a tank-to-jet diameter ratio of 25, 
and Sterbentz’s (ref. 1 1) used a ratio of 21. In practice, the jet diameter could be 
optimized by trading the kinetic energy reduction associated with a l^ger-diMiieter 
nozzle a gain st the mas s and volume of the device, but a diameter ratio of 25 is expected 
to be representative of the optimum range for most cases. 

The fill firaction chosen for the experiment is approximately 85%, representing a nearly 
full propellant tank (tihe actual fill fraction as flown was 83%). Ideally a variable fill 
fraction would be tested, but for simplicity, reliability, safety, and practicality a smgle 
fill level had to be selected, and 85% was chosen to represent a nearly-fuU condition. A 
high fill level (75% to 95%) may be a critical design condition for pressure control 
systems for several reasons: (1) the rate of pressure rise is greatest with a small ullage 
volume, for a given rate of heat input to the ullage and/or rate of evaporation; (2) the 
greater the fill fraction is, the further the liquid jet noay have to travel to reach the ullage, 
therefore the greater the required flow rate to produce a given jet Weber number (or 
noixing effectiveness); and (3) for a non-vent^ fill of a tank in low-g, the nearly-full 
condition is probably the fill rate-limiting condition because the interfacial areas available 
for ullage cooling and condensation are small, hence the effectiveness of the mixing 
could strongly control the fill time and/or the final achievable fill level. 

The experiment operation consists of a series of test rans in which the pressure is raised 
by addition of a fixed amount of heat and then is reduced by mixing. This intermittent 
mixer operation is to simulate an on-demand pressure control system, which could be 
activated whenever the pressure in a storage tank rises to a predetermined level and 
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operated until the pressure falls back below the fixed control band. However, TPCE 
uses a fixed duration for each heating and mixing phase, for simplicity and reliability, 
rather than operating "on-demand". Each run consists of a ten-minute pressure rise 
period (unless the pressure increases by 21 kPa [3 psi], in which case the heating phase 
is terminated early for safety reasons). This is followed by a fifteen-minute pressure 
reduction (mixing) phase. Another fifteen minutes with no activity is added to each run 
to allow fluid motions to damp out, and the equilibrium liquidA^apor interface shape to 
be reestablished, before beginning the next run. 

Heat leaks into the fluid, and the resulting temperature gradients, are simulated by two 
heaters immersed in the fluid. The heating rate is high enough to create a significant 
pressure rise within a short time for each test run, but low enough (<25 watts) to avoid 
potential in-flight hazards associated with overheating the refrigerant vapor. To help 
meet safety requirements the heaters are designed to have sufficient thermal nmss to 
avoid the possibility of their reaching the decomposition ten^rature in ten minutes of 
operation at 25W. A surface area of approximately 130 cm2^(20 in2) is chosen (see the 
Experiment Apparatus and Operation section for the actual values) to provide a heat flux 
of up to 0.2 W/cm2 (1.25 W/in^). This heat flux is below the expected burnout heat 
flux of 0.5 W/cm2 (3.5 Wfin^) at 10-^g (calculated firom an equation by Zuber, ref. 19), 
so that film boiling is not expected to be encountered (film boiling would reduce the 
amount of heat added to the fluid due to its lower heat transfer coefficients, so nucleate 
boiling was desired). At a 25W heating rate, a maximum pressure rise of 14 to 21 kPa 
(2 to 3 psi) over a ten-minute heating period in low g was predicted by a lumped 
thermodynamic model of the fluid that accounts for evaporation, condensation, and heat 
conduction. 

Each heater is near a tank wall but offset at least 0.6 cm (0.25 in) so that the majority of 
the heat flows into the fluid rather than out of the tank. One heater, designated heater A, 
is at the end of the tank opposite the mixer nozzle. The second heater, heater B, is 
located off the tank major axis along the side of the tank. Heater B's location is intended 
to produce an asymmetric initial condition (i.e., the bubble(s) and the temperature 
gradient are not centered on the path of the mixing jet) so that the effects of a realistic 
asymmetric, three-dimensional initial condition can be assessed. Heater B is also offset 
a greater distance from the tank wall to ensure that it, and any boiling activity, is visible 
in the photography. 

A range of mixer flow rates is selected that spans a ten-to-one range, to include both 
penetrating and nonpenetrating flow regimes. The flow pattern conelation of Aydelott 
(ref. 15) was used in selecting the values in the test matrix, given in figure 4. The figure 
shows the flow rates chosen along with the corresponding flow patterns predicted by 
Aydelott's correlation. Four "primary" flow rates are the highest priority: two Pattern 
IV points, one Pattern I point, and zero flow (to assess self-mixing, or the pressure 
decay rate without forced circulation). Each primary test point is run three times with 
each of the two heaters to assess the repeatibility of the results, and is also run once with 
both heaters on (with the heaters connected in series so that the total heat input is one- 
half the normal rate, and one-quarter the normal heat flux). Five "secondary" flow rates 
are chosen to extend the range above and below the primary points and to test 
intermediate values for which Patterns II and HI are predicted. Each of these is mn once 
with each heater. There are thus 38 total test runs: 24 primary arid 14 seconda^. The 
runs are arranged in a pre-planned randomized order, with the primary runs being 
executed first. Each run lasts for approximately 40 minutes, for a total test duration of 
25.3 hours. 
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Figure 4: Test Matrix and Predicted Flow Patterns 

The acceleration requirements for the experiment operation are consistent with normal 
STS orbital operations. Steady-state g-levels of l(Hg ot less are sufficient to produce a 
capillary force dominated fluid raientation, and are low enough to have a negligible 
effect on the jet geyser. This g-level also is low enough that buoyancy-driven free 
convection is negligible. It is desirable that for some of the test runs Ae liquid be settled 
at the mixer end of the tank and the ullage bubble at the opposite end — this axisymmetric 
fluid configuration is most similar to i^vious experimentatitm, most amenable to 
analysis, and of interest for those applications in which orbital drag will keep the liquid 
settled at one end of the tank. This initial condition is provided by cnienting the 
experiment tank with the nozzle end aft, such that the launch and OMS insertion bums 
will position the liquid at that end. Furthermore, a period of "tail-first” Qrbiter 
orientation is specified so that on-orbit atmospheric drag tends to settle the liquid at the 
aft end (a period of eight or more hours of tail-first time was requested; nearly 12 hours 
was actually provided on STS-43 as will be discussed later). In addition to this tail-first 
attitude, the more typical "nose-first" orientation and other attitudes were also desired for 
portions of the test timeline to assess the effect of initial fluid orientation on the results. 

Orbiter reaction control system (RCS) thruster firings, crew activities, and other events 
can impart brief acceleration levels of up to lO'^g. TOs noise, above and beyond the 
steady-state acceleration due to drag, can induce significant transient motions in the fluid 
and the interface. If their duration is sufficient, these events can cause the fluid to be 
reodented. The level of noise has been shown in previous flights to be greatly 
dependent on the phase of the mission and the activities of the crew and other payloads. 
O^ration of the experiment over at least a 24-hour period provittes some data during 
periods of high activity and some during a sleeping period. This provides an indication 
of the effect of acceleration disturbances on pressure control system performance, 
including self-mixing. 

Instrumentation 

The primary data to be obtained are visual: the fluid orientations and the jet flow 
patterns as functions of flow rate. These are recorded by video cameras that photograph 
the fluid through transparent tank walls. In each run the first two minutes of the heating 
stage are photographed, so that the fluid initial conditions and the onset of boiling can 
observed. The camera is then turned off until the beginning of the mixing stage, when it 
is turned on again for four minutes of photography. There is thaefore a total of nearly 
four hours of real-time video (the maximum capacity of the tapes used). As aids in 
correlating the video data with the experiment timeline and digital data, brief periods of 
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camera operation without lamps provide separation of the runs and of the heating and 
mixing stages, and two light-emitting diodes (LEDs) flash a code at the beginning of 
each run that identifies the run number. The LEDs also flash once per second during 
photography to aid in timing any events seen on the videotapes. 

Quantitative data taken are the tank pressure, the mixer flow rate, the temperature at four 
locations in the fluid, the temperature of each heater plate, the temperature of the fluid in 
the mixer flow, the heater power, and the acceleration in each of tluee axes. The 
purpose of measuring accelerations is to characterize the relative noisiness of different 
portions of the test, and to record the timing, magnitude, and duration of acceleration 
events that may influence the experimental results and possibly explain behavior seen in 
the video or pressure data. Determination of the precise value of the orbital drag 
acceleration is not a requirement of the experiment because, except fOT its effect on the 
initial ullage location, the steady-state g-level is not high enough to affect mixing 
phenomena (whereas the acceleration spikes could). 

Specifics on the instrumentation (range, resolution, and accuracy) are given in the 
following sections. 
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EXPERIMENT APPARATUS AND OPERATION 


STS Accommodations 

The Tank Pressure Control Experiment is installed aboard the Space Transportation 
System (STS) in a Get-Away Special (GAS) carrier system. The GAS (or Small Self- 
Contained Payload) earner provides standardized intrafaces with the STS and can 
facilitate flight manifesting, integration, and safety. Payloads using the GAS carrier 
must be self-contained and autonomous, requiring no (Arbiter power or data services and 
minimal crew control. 

TPCE weighs 84 kg (186 Ibm) and fills the 0.14 m^ (5 ft^) GAS canister. For the STS- 
43 mission the experiment was mounted on the port side of the payload bay, on the sill, 
approximately one-third of the way aft of the forward end of die bay. TPC^'s GAS 
canister is purged and sealed, and contains one atmosphere of dry nitrogen gas 
throughout the flight The exp«iment is controlled by a microcomputer that activates 
pumps, heaters, lamps, and instrumentation in a pre-programmed sequence. The 
computer also record data and stores it in solid-state memory. All payload power is 
from internal alkaline batteries. 

Figi^ 5 is a photograph of the payload being installed in a GAS canister. Figure 6 is a 
sectional drawing of TPCE within the canister showing major elements of the apparatus 
(most details are omitted or simplified in this drawing for clarity). 

Operation of the payload is initiated by a baroswitch, which is GAS-standard 
equipment The baroswitch senses the depressurization of the payload bay during 
laimch and closes a relay which connects the payload's battery pack. Operation of the 
experiment is autonomous after that point As a backup to the baroswitch, the STS-43 
crew checked the status of the relay to ensure that it was closed. After allowing 
sufficient time for the experiment to complete its planned sequence, the crew opened the 
relay via the autonomous payload controller (APQ in the aft flight deck, removing 
power from the payload. The baroswitch, power control relays, and ATC interfaces are 
contained in the GAS electronics area visible in figiue 6. 

A nnoderate thermal attitude was provided during the operating time of TPCE. This 
included a predominantly bay-toward-Earth orientation, with brief excursions toward the 
sun and deep space. The GAS canister is well insulate^ and the payload itself is 
thermally isolated from the canister, so the thermal environment for the experiment was 
stable and benign. A "tail-first" Orbiter attitude was requested for at least eight hours 
during experiment operation so that orbital atmospheric drag would tend to setde the 
liquid at the aft end of the tank. Nearly twelve hours of this attitude were actually 
provided, which included the first sleeping period. 
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Figure 5 ; TPCE and Get-Away Special Carrier 
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0.826 m (32.5 in) 
container length 






structure 


All components of the payload arc mounted inside a composite shell. The composite, a 
foam-fiberglass-epoxy sandwich, provides high strength and stiffness plus low thermal 
conductance to the environment. It is lightweight and low-cost, and allows flexibility 
for mounting components during assembly. The shell has an octagonal cross section 
and is constmcted in two pieces: a main "lower" section, and an "upper" tray. The two 
sections arc mounted to each other and to the GAS canister's upper plate by eight tie 
bars, as shown in figtue 7. Four equally- spaced adjustable bumpers carry lateral loads 
to the inside of the GAS canister at the lower end. 

The fluid tank, the heaviest assembly in the payload, is mounted to six aluminum 
brackets via neoprcne grommets. The grommets provide further isolation of the tank 
from heat conduction, vibration, and stresses. The aluminum brackets arc attached to 
the lower section of the composite shell. Figure 8 is a photograph of the lower end of 
the payload that shows some of these details. For this photograph the end panel was 
removed — normally this panel closes off the bottom of the payload and provides a 
reflective interior surface for lighting. Figure 9 is a drawing showing the installation of 
the tank as viewed from the "upper" end of the payload (tfie camera’s viewpoint). It 
also shows the locations of the bumpers and of some of the fluid system components. 

The upper composite tray contains the battery pack, computer, and video cameras. It is 
removable for access to the fluid system, memory modules, and electronics contained in 
the lower section. 














Figure 9: Tank Mounting Orientation 


Fluid System 

In accordance with the technical requirements stated previously, the tank is constructed 
in the approximate shape of a 25.4-cm (10-inch) diameter cylinder with hemispherical 
ends and a total inside length of 35.6 cm (14 inches). Its total calculated volume is 13.7 
liters (0.485 ft^). Its capacity was measured to be 21.8 +/-0.5 kg(48 4 -/- 1.0 Ibm) of 
refrigerant 113, or 13.9 liters (0.49 ft^). This measured capacity includes the volume of 
fluid contained in the pumped lot^ but not the small volume of the pressure transducer 
manifold. 

The tank itself is constructed in three main pieces: two transparent domes and a stainless 
steel oval-shaped center ring. The two domes are thennoformed ftom 0.1 -cm (0.25 -in) 
original thickness sheets of a crazing-resistant grade of acrylic plastic. They are bolted 
to the steel ring and are sealed against it with buna-N O-rings. The gaps between the 
flange bend radius and the flat surfaces of the steel ring are filled with epoxy, as shown 
in figure 10. The epoxy filler provides a smooth, rounded sinface to avoid flow 
disturbances. All attachments and penetrations to the tank are made through the steel 
ring, including plumbing, thermistor probes, heater supports and feedthroughs, and 
LAD mounts. Figure 1 1 is a photograph of the tank assembly, and figure 12 shows its 
installation in the payload from the point of view of the cameras. All three fluid tubes 
leading to the tank are flexible metal hose which, in addition to the neoprene tpk 
mounting grommets, provides strain relief and reduces transmission of vibrations to the 
fluid. 
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Figure 11: Tank Assembly 
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Figure 12: Fluid System Installation 

Figure 13 is a schematic of the fluid system. The tank contains six temperature probes 
consisting of thermistors encapsulated in stainless steel sheaths. Thermistors T2, T4, 
and T5 reach 2.54 +/-0. 15 cm (1.0 +/-0.06 in) from the inside wall of the tank, and 
thermistor T1 extends 1.27 +/-0.15 cm (0.5 4-/-0.06 in). Thermistors T3 and T6 are 
attached to heaters A and B. Thermistor T7 is another probe inserted into the mixer flow 
downstream from the flow meter. A pressure transducer is connected to the tank via a 
small manifold and tube. 

Also connected to the manifold are a pressure relief valve and a pressure switch. A 
second relief valve is mounted directly to the side of the tank. The computer is 
programmed to turn the heaters off if the pressure reaches 157 kPa absolute (22.8 psia), 
which in the nominally 103-kPa (15-psi) GAS canister environment equates to 54 kPa 
gage (7.8 psig). As a backup, the pressme switch is a safety device to interrupt power 
to the heaters in the case the pressure rises above 62 kPa gage (9.0 psig). In addition, 
the relief valves prevent any potential for tank rupture by venting fluid if the pressure 
reaches 69 kPa gage (10.0 psig). The tank's maximum design pressure is 86.1 kPa 
(12.5 psig), with a design ultimate safety factor of 2.5. (Under normal operating 
conditions the tank was not expected to reach a positive differential pressure at any time 
during the STS-43 mission, and this proved to be the case.) 

Two pumps plumbed in parallel are used to provide the mixing flow. A single pump is 
used at lower flow rates (with the capability of switching to the other pump in case the 
first pump fails to deliver flow). Both are operated in the runs requiring higher flow 
rates. To reduce the net positive suction head requirements, the pump suction line and 
inlet ports are designed for minimum pressure loss. 
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The pumps arc supplied with vapor-frce liquid at all times by the LAD, which is a 
semicircular channel made of stainless steel plate. The side of this channel facing the 
tank wall is micro-etched with 0.38-mm (0.015-in) diameter perforations and has a 21% 
open area fraction. The gap width between this porous surface and the tank wall is 0.6 
cm (1/4 in). 

Two heaters are immCTsed in the fluid. Heater A is located 0.5 cm (0.2 in) from one end 
of the tank (at its maximum distance), and heater B is located off the tank centerline 
approximately 2.5 cm (1.0 in) from the side wall. Heater A is bent to a 12.1-cm (4.75- 
in) radius to follow the curvature of the tank wall, while heater B is flat The two 
heaters are constructed of an etched-foil element encased in silicone rubber insulation, 
which is sandwiched between two 304L stainless steel plates as shown in figure 14. 

The perimeter of the sandwich is welded to seal the element from the fluid. The heating 
element is a 9.6 cm by 6.6 cm (3.8 by 2.6 in) rectangle with a 0.64-cm (0.25-in) 
circular cutout so that thermistor probe T4 can pass l^ugh the heater. Each steel plate 
has a thickness of 0.191 cm (0.075 in), narrowing to 0.089 cm (0.035 in) near the 
perimeter weld. The outside dimensions of the assembly are 10.46 by 7.42 cm (4. 12 by 
2.92 in), for an overall surface area (both sides included) of 155 cm^ (24.1 in^). This 
area includes the small perimeter flange that extends 0.30 to 0.41 cm (0.12 to 0.16 in) 
beyond the element. The heater assembly is welded to a standoff tube which supports 
the assembly and contains the leads. Ad^tional support is provided by a thermistor 
probe bond^ (using thermally-conductive epoxy) onto a small area of greater thickness. 
The total mass of each heater, excluding the standoff tube and thermistor, is 0.214 kg 
(0.472 Ibm), and the thermal capacitance is estimated to be 0.10 kJ/K. Figure 15 is a 
photograph of the two heaters. 
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Figure 14: Heater Construction (not to scale) 








Power Supply 

The payload is powered by an internal battery pack consisting of 96 F-size alkaline 
(non-rechargeable) cells. The cells are connected in a series-parallel arrangement that 
provides 12 and 18 volt (nominal) outputs. The total battery pack capacity is 20(X) watt- 
hours under constant drain rates, which is approximately three times the energy required 
by the payload under normal conditions. The extra capacity provides voltage stabdity 
under peak drain rates, better low-tempea-ature performance, and reserve capacity in case 
of partial battery failure or other paylc^ failures. The cells are pac^ged into 48 two- 
cell stacks which are glued into a rigid polyurethane foam block imide a fiberglass shell. 
This forms a mechanically strong block with minimum mtemal void volume. A sn^ 
volume at the top of the pack, containing the wiring, fuses, and protective diodes, is 
packed with absorbent insulating felt material and covered with a fiberglass lid. The 
pack is designed to contain any liquid electrolyte that might leak in the event of damage, 
but without containing any gases that might bie evolved. At 22 kg (48 Ibm), the battery 
pack is a large portion of the payload mass. Figure 16 shows the battery pack installed 
in the upper tray along with the computer and camera enclosure. 
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Figure 16: View of Payload Showing Battery Pack, Computer, Cameras 
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The ground terminal of the battery pack is coimected to the GAS container, which is 
grounded to the Orbiter. The 12V and 18V power leads are connected to the payload via 
the GAS power relays so that all payload power is controlled by these relays. Payload 
lamps and cooling fans are supplied with unregulated 12V DC power, and the pumps 
and heaters use unregulated 18 V DC. The mixing punqjs are controlled by a pulse- 
modulation circuit that also compensates for the d^lining voltage as the batteries 
discharge, to produce the desired pump speed and flow rate for each test run. Voltage 
regulators also supply regulated power at +15, -15, +5, +6, and +8 volts to the 
computer, instrumentation, and video cameras. The power regulating circuits and other 
components that dissipate significant amounts of heat are located on a circuit board in the 
main compartment of the payload. Two small cooling fans circulate the nitrogen purge 
gas in the payload to provide forced-convective cooluig of these components. 

Controls and Data Storage 

A microcomputer mounted in the upper tray controls the payload operation and records 
the output of the instraments. The computer, visible in figure 16, is based on the ultra- 
low-power C-44 bus. It uses mostly off-the-shelf complementary metal oxide 
semiconductor (CMOS) boards having low energy consumption (and consequently low 
heat dissipation). Four of these circuit boards are housed in a card cage: (1) the central 
processing unit (CPU) board which uses the 80C85 processor, (2) a 16-channel analog- 
to-digital ccmverter board, (3) a custom-built output control and equipment interface 
board, and (4) a memory interface and control board. Figtue 17 is a block diagram of 
the control system, A custom-built circuit board, called Ae conditioner board, contains 
the interfaces between the computer and the instrumentation (amplifiers, voltage 
dividers, pulse counter, filters, and other signal conditioning circuitry) as well as some 
of the output interfaces and voltage regulators. 



Figure 1 7: Controls and Data Storage Block Diagram 

Software for the controller is written in 8085 assembly language. The code is stored in 
a non-erasable, fusible-link type of read-only memory so that the instmctions cannot be 
altered by cosmic radiation. Since CMOS devices are known to be susceptible to upset 
by the taxation environment within a GAS canister in flight, critical variables that must 
be stored in random access memory during operation are stored in redundant locations. 


33 




This makes it possible to recover proper data values (such as the program counter) in the 
event that a value is erased or altered by radiation. 

The memory interface board has its own CPU to ccMitrol the storage of digital data. It is 
connected to three solid-state memory nKxlules, each of which stores (me megabyte of 
data on nonvolatile erasable programmable read-only memory (EPROM) chips. These 
three memory modules are mounted in the payload's lower section and are removable 
for extraction of the data. Nearly three megabytes of data are stored during execution of 
a normal mission. Following the mission or a mission simulation test the data is 
downloaded to a mainframe computer for reduction into engineering units, application of 
calibration data, and plotting. 

Instrumentation 

The primary experiment instruments include a flow meter downstream of the pumps, a 
pressure transducer, the seven fluid temperature sensors (thermistors) discussed 
previously, and three accelerometers. Figure 18 summarizes the range and accuracy of 
each instrument, plus the rate at which dato are taken during each stage of a test run. 

The atxairacy figures are estimated based (mi the manufacturers' stated accuracies and on 
their performance in flight and during pre- and post-flight tests. The pressure 
transducer, which measures absolute pressure, was found in post-flight testing to have a 
shifted bias value, i.e. all readings were offset by a snoall constant vdue. All pressure 
data were subsequently corrected for this bias shift, which was likely caused by torque 
on the transducer fitting during installatitm. Each "thermistor" is actually a set of th^ 
matched thermistors, packaged together in the tip of the probe, that are read in such a 
way as to have a linear output Self-heating errors are muiimized by applying current to 
each thermistor only during the brief period that a reading is taken. 


INSTRUMENT 

RANGE 

RESO- 

LUTION 

ACCURACY 

SAN 

IPLERATE 1 

Stage 

2 

stage 

3 

stage 

4 

pressure transducer 
(absolute) 

0tol73kPa 
(0 to 25 psia) 

0.04 kPa 
(0.006 psi) 

0.35 kPa 
(0.05 psIa) 

3/sec 

3/sec 

1/sec 

thermistors (7) 

-20 to +70‘C 
(-4 to +158‘F) 

0.02*C 

(0.036'F) 

0.1 *c 

(0.1 8*F) 

1/sec 

1/sec 

1/min 

flow meter 

0 to 3.8 l/min 
(0 to1.0gpm) 

0.02 l/min 
(0.005 gpm) 

0.04 l/min 
(0.01 gpm) 

- 

1/sec 

- 

heater power 

0-30 W 

0.1 W 

0,1 w 

1/sec 

- 

- 

accelerometers (3) 

+/- 0.01 g 

2.4xl0"*g 

5 X 10*^g abs. 
1.2x10^g rol. 

3/sec 

3/sec 

3/sec 


Figure 18: Instrumentation Summary 

The accuracy of the three accelerometers is given for absolute readings and relative 
readings. The bias (zero-point) reading of these units is subject to drift over time and to 
temperature effects, altering the absolute readings. The scale factor, which governs 
relative readings such as the magnitude of the acceleration produced by a thmster firing 
versus the background level, is less susceptible to these errors. The X-axis 
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accelerometer, being of highest interest for axial jet mixing, is thermally modelled so that 
its output is adjusted in the post-flight data reduction for changes in both its bias value 
and scale factor as functions of temperature. The ouqjuts of all three accelerometers are 
filtered, prior to sampling by the computer's analog-to-digital converters, to remove 
high-ftequency components of the signal. The an^og filter used has a time constant of 
0.375 seconds. 


Other data recorded, for diagnostics and to perform the temperature compensation of the 
X-accelerometer, are two additional temperatures: (1) inside the X-axis accelerometer, 
and (2) on the analog-to-digital circuit board in the main computer housing. These have 
accuracies of approximately 2*C (3.6“F) and are recorded once per run. Voltage-divider 
circuits are also used to sense the heater power and the battery voltages. The voltage at 
the 18V nominal battery terminal is recorded once per minute with an accuracy of 
approximately 0.18 V. The 12V level is also measured once per minute to determine that 
the cells are not overdischarging, but this data is not recorded 

The three accelerometers are installed on one of the tank supports such that their axes are 
aligned with the STS Orbiter coordinates. Figure 19 shows die TPCE acceleration axes 
witfi respect to the fluid tank, to the Orbiter, and to the picture visible when viewing the 
videotapes. 


TPCE Coordinates: 



Relative to Video Image: 


(■fZaway 


+Y 



from viewer)! 
+X 



Figure 19: Accelerometer Axes 


Photography 

Video photography and recording is accomplished with two modified commercial 
camer^recorder devices. These compact consumer-type 8mm format camcorders 
provide two hours of photography each. They are mo^ed to use payload power 
(supplied at 8.0 volts) with computer control of their functions. Certain portions of the 
cameras are removed and severe switches are locked into position. They are mounted 
in a stainless-steel enclosure with approximately 0.6 cm (0.25 in) of foam padding on all 
sides and are allowed to "float" inside it for their protection from vibration and shock 
during launch. Wide-angle adaptor lenses are mounted on the enclosure, which is 
airtight to eliminate concerns atout outgassing of lubricants and other materials in the 
camcorders. Figure 20 is a photograph of the enclosure. 

The cameras view the tank through an opening in the upper tray. Their fields of view 
are similar, with lines of sight paj^el to the centerline of the payload, except that the 
right-hand camera image is shifted slightly towards heater B to ensure that boiling 
phenomena on the heater plate are visible. 

Four 12-watt halogen incandescant lamps provide illumination of the tank. Their light is 
directed at the panel at the lower end of Ae payload. This panel is reflective (a glass- 
beaded aluminum finish) with a 0.64-cm (0.25-in) square black grid painted on it. The 
grid's purpose is to enhance visualization of the liquid/vapor interface — ^the difference 
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Figure 20: Video Camera/Recorder Assembly 

between the refractive indices of the phases produces discontmuities in the grid lines at 
the liquid/vapor interfaces. The fluid is tint^ with a trace amount (<0.5 ppm) of orange 
dye (Sudan IV) to further aid in discriiriinating liquid from vapor in the picture. 

The lamps are sealed with glass plates (which are frosted to diffuse the light), preventing 
contact between the hot bulbs and any leaked refrigerant, which could create a 
decomposition hazard. Each bulb assembly is mounted to an aluminum bracket which 
blocks most direct viewing of the tank by die bulb and which also acts as a heat sink to 
minimiyf* radiative heating of the fluid. Each bracket is mounted to a tank supper^ 
providing a further heat sink, and the negative electrical lead of each bulb (which is 
connected to the hot filament) is thermally shorted to the heat sink using a large-diameter 
ground wire. 
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Operating Timeline 

Depressurization of the Qrbiter payload bay during ascent causes the baroswitch to 
activate the payload approximately one minute aftCT lift-off. At this time the computer 
be^s to operate and starts a six-hour delay period. This delay allows orbit to be 
achieved and the payload bay doors to be opened before starting the experiment After 
this delay period the computer turns on power to the instruments, allows a four-minute 
warmup period, and performs a ten-second startup health test. After checking for low 
fluid pressure, which would cause execution of a fluid warmup contingency operation, 
and for high fluid pressure, which would trigger a high-pressure contingency, the test 
matrix of thirty-eight mns is begun. 

Figxu« 21 shows the order of the tests and the camera, pump, and heater assigned to 
each. A run begins by first checking for excessively lugh pressiue and, if found, 
attempting to reduce it befOTe turning on the heater(s). "niis is identified as Stage 1 and 
its duration is one-third second unless the high pressure contingency is encountered. 
Stage 2 is the heating stage and its total duration is ten minutes, unless the pressure rises 
more than 21 kPa (3.0 psi) before that time. Data-taking begins at tfie start of this stage, 
and video photo^phy begins at approximately six seconds. First a code is flashed on 
the two L^s to identify Ae run number, then the lamps are turned on at 13.33 seconds 
from start of the run. One second later the heater is turned on. The lamps and camera 
are turned off at two minutes, and at 10:(X).33 the heater is turned off. 

Stage 3 is the mixing stage. At 10:02.33 the camera is again turned on and, after 
allowing a brief period of black picture (to distinguish the stages in the video), the lamps 
are turned on at 10:06.33. Two seconds later the pump(s) is turned on. The camera is 
turned off after four minutes, and the pump(s) remains on until 25:00.33, the end of 
Stage 3. Stage 4, the recovery phase, runs for fifteen minutes with data taken at a 
reduced frequency (except for accelerometCT data, which is recorded at 3 Hz 
continuously). No heating, mixing, or photography occurs during Stage 4. The entire 
duration of each mn is 40:01, and the next run begins immediately afterward. 

The complete sequence of test runs lasts over 25 hours (25:20:38), not including the six- 
hour initi^ delay and warmup period. The test matrix concludes at approximately 3 1 .5 
hours after lift-off. At this time the computer removes power finom most of the payload 
con^nents and puts itself into a low-power sleep mode. The STS-43 crew turned off 
the power relay two days after launch, which allowed sufficient time for this nominal 
sequence plus all of the programmed contingencies. 
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LOW-G AIRCRAFT TEST 


Aircraft Test Objectives and Approach 

A prototype payload was constructed for a low-gravity test on the NASA Lewis 
Research Center Leaijet prior to assembly of the STS payload. Up to 30 seconds at a 
time of reduced gravity are obtained during parabolic ^ghts of this aircraft. The 
purpose of this test was to reduce technical risk for the STS test by performing a 
comprehensive systems checkout in a low-g flight environment. 

The general objective of the aircraft flight test was to verify the overall operation of the 
payload inside a GAS canister in a low-g flight environment It was intended to uncover 
any possible unexpected performance problems that are not seen in normal gravity 
testing. In addition, the test would provide experience in shipping the payload, 
integrating it with the GAS canister, deintegrating, and extracting data after the flight. 
Five specific objectives were identified for the aircraft flight test: 

1 . Flow Patterns. The flow rates for TPCE were selected using a correlation that is 
based on subscale drop tower tests. Confirmation that these flow rates will span 
the desired range of flow patterns was desirable to maximize the value of the data 
returned from the STS flight. 

2. Cavitation. Proper operation of the pumps in the absence of a gravitational head 
pressure was a concern. The low-g test would verify that punq) cavitation is 
suppressed when 7-14 kPa (1-2 psi) of noncondensible gas is added to the 
ullage to provide net positive suction pressure (NPSP). 

3. Fluid Foaming. Another concern was the possibility that the jet mixer would 
cause excessive bubble formation, particularly at the highest flow rates, resulting 
in foaming and frothing of the fluid in the tank. If the fioth were to persist, 
rather than coalesce, subsequent test runs would be compromised. There was 
insufficient experience with bubble coalescence and collapse times in low g to 
completely exclude this possibility without a test 

4. Heater Operation. Observation of the heaters was desired to possibly determine 
whether nucleate boiling is obtained at their heat flux in low g. This could 
indicate whether the heaters are capable of providing the desired pressure rise 
during the heating phase of the experiment (It was recognized however that the 
short duration of low-g time available might not permit boiling to be established 
in the aircraft test) 

5 . Photographic Coverage. Because the liquid and the tank walls reflect light the 
field of view and the ability to see detail in the tank are highly affected by the 
liquid/vapor orientation. Since a low-g fluid orientation, particularly a curved 
interface, is difficult to achieve in the laboratory, it was necessary to test the 
adequacy of the video photography optical design and lighting on the aircraft. 

The prototype payload was installed inside a GAS canister and hard-mounted to the 
cabin floor of Ae Leaijet. The payload was controlled in real time by on-board 
experimenters via a small laptop control box. A simplified version of the flight software 
was installed in the payload computer which caused it to record data at a constant rate 
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while permitting manual control of the pumps, heaters, lamps, and cameras. The 
payload was powered by its internal battery pack, and it recorded data into its internal 
memory for post-flight extraction, although several channels were monitored during the 
flight via test ports on the control box. Both video cameras operated continuously 
during each flight test The payload was oriented so that the tank's longitudinal axis 
was vertical during level flight, with the mixer nozzle and LAD at the bottom. This 
ensured that the mixer pump was supplied with liquid during the parabola entry 
maneuver and that the LAD was initi^y filled with liquid. 

Twenty-four parabolas were flown during four flights. During the first 12 runs the tank 
was pressurized with helium gas to a tot^ NPSP of 15-21 kPa (2.2-3 psi). For the last 
12 runs the NPSP was reduced to approximately 6 kPa (1 psi). Each of the flow rates, 
including primary and secondary points, from the STS test matrix was tested several 
times. In each run the heater was turned on during level flight, and the pump(s) was 
adjusted to produce the desired indicated flow rate and then turned off. Sever^ seconds 
after entering the low-g portion of the maneuver the pump(s) was turned on and was left 
on until the 2.5g pull-up began. In the last three runs the pump(s) was turned on prior 
to low-g entry and left on throughout the parabola. 

Aircraft Test Results and Conclusions 

The payload operated well throughout all four flights. All fli^t data were extracted and 
reviewed. It was discovered that in the first flight (the first six test runs) the pilot’s 
display of the target acceleration, which is used for pitch control during the parabola, 
was inadvertently set for 0.05g rather than Og. This resulted in a flat liquid/vapor 
interface with no observed jet penetration during those six runs. The remaining flights 
provided acceleration levels which were approximately centered around Og in all three 
axes. However there were sigitificant excursions from Og which frequentiy reached 
full-scale on the accelerometers (+/-0.01g) in all runs. Figure 22 shows typical 
accelerometer readings from the three channels. (In this plot, positive "X" is toward the 
starboard wing, or away from heater B; positive "Y" is toward the nose of the aircraft, 
or toward the cameras; and positive "Z" is "up" during level flight, parallel to the 
direction of jet flow out of the nozzle.) 

In reviewing the video data it was evident that aircraft tests cannot substitute for in-space 
operation of an experiment of this kind. The quality oi the acceleration environment was 
inadequate to prxxiuce a surface-tension-dominated fluid orientation, and there was 
insufficient time in each parabola for the liquid motions induced by the pull-up (and 
amplified by the g-level reduction) to damp out Significant velocities remained in the 
liquid in each run and made it difficult to establish a mixer flow pattern. However, it 
was determined that the highest flow rates occasionally appeared to penetrate the ullage, 
and the lowest flow rates dad not appear to be capable of j^netrating the surface. Since 
the planned test matrix therefore appeared to span the desired range of flow patterns, it 
was concluded that the matrix was validated. 

Boiling was not observed at the heaters in any run. With the significant residual fluid 
velocities present during the maneuvers, the low heat flux was insufficient to overcome 
the forced-convection heat transfer at the heater plates. Consequently specific objective 
number 4 was determined to be unobtainable with this kind of test 
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Figure 22 : Aircrcfi Test Accelerometer Traces: One Typical Parabola 

The other objectives were successfully met It was confirmed that cavitation would be 
suppressed by the level of noncondensible gas plaimed to be added to the tank. Some 
cavitation — ^manifested as a reduction in flow rate during low-g — ^was seen at the lower 
value of NPSP, but the experiment was otherwise unaffected. Foaming of the liquid 
was also determined not to be a valid concern. Even in this worst-case tes^ ^th 
considerable sloshing of the fluid in combination with the high flow me mixing, the 
bubbles were seen to rapidly coalesce and no tendency to form a persistent froth was 
found. The photographic quality proved to be goo<k providing visualization sufficient 
to meet the experiment requirements under all conditions. 

The experiences resulted in several improvements to the payload design and operations 
that proved valuable on the STS flight One camera was re-aimed to yield better 
coverage of boiling at heater B. The burners were redesigned, and a low-rate leakage 
problem with the tank assembly was identified and solved. Better procedures for 
erasing and downloading from the memory modules were develop^ to reduce the risk 
of data loss. No other problems relating to operation of the apparatus in low gravity 
were found. 
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IN-SPACE TEST 


Payload Performance 

General. The TTCE payload was flown in space on Shuttle mission STS-43 from 
August 2 to 1 1, 1991. Eight hours after launch , the Shutde crew verified that the 
baroswitch had activated the payload. The "tail-first" (-Xw) Orbiter attitude was 
initiated at approximately day zero, seven hours and fifty minutes (0/07:50:00) mission- 
elapsed time (MET), and continued through approximately 0/22:0OKX) MET, for a total 
duration of alx)ut 14 hours. TPCE was deactivated on schedule two days after launch at 
2A)1:32:48 MET. Thus TPCE was powered up for approximately 49.5 hours. 

After the flight the payload timeline was correlated with the mission timeline by 
comparison of mission events recorded on the accelerometer data channels. The 
deployment of the primary payload, a Tracking and Data Relay Satellite and its Inertial 
Upper Stage (TDRS/IUS), is visible as a 3-second sawtooth- wave centered at 0:07:09 
on TPCE’s clock. The actual time of deployment of the TDRS/IUS payload is known to 
be 6:12:39 MET, so that TPCE's time-zero is at 06:05:30 MET. 

After recovery of the payload it was inspected and tested. No damage was observed to 
any components and no failures occurred. The systems testing performed at Kennedy 
Space Center showed no anomalies. The liquid level was inspected visually relative to 
references on the tank and no change was observed. After the ullage gas was sampled 
for analysis, the tank was drained and weighed. Allowing for evaporation during the 
sampling, the original fill quantity of 18.1 +/-0.5 kg (40 +/-1 Ibm) was confirmed. This 
corresponds to a fill fraction (at room temperature) of 83%. 

Ullage Gas Makeup. Measurements prior to shipment showed that pressure 
changes were occurring over time in the fluid system due to the permeability of the 
polj^eric tank and seds. It was originally intended that the ullage consist only of 
refingerant-1 13 vapor plus 14 kPa (2.0 psi) partial pressure of helium. However, the 
total partial pressure of noncondensible gases (Pnc) was observed to be dropping slowly 
during storage prior to the flight Pnc is defined as the difference between the total 
pressiue (with the tank well-mixed) and the vapor pressure of refrigerant 1 13 as 
calculated for its temperature at the time of the measurement Based on the flight data 
(i.e., the pressures and temperatures recorded when the fluid was well mixed), the value 
of Pnc at the time of experiment operation was approximately 4.1 +/-1.0 kPa (0.60 4-/- 
0.14 psi). 

With the reduced Pnc to provide pump suction pressure (4.1 kPa versus the planned 14) 
there was concern that cavitation woitid limit the pumps' ability to provide the planned 
flow rates. In flight, however, for the two runs having the highest flow rate (runs 37 
and 38) the measiued flow rates averaged nearly 90% of the planned values. The other 
two effects of cavitation that had been observed in ground testing, besides flow rate 
limiting, were bubbles in the outflow line and erratic flow meter data due to bubbles in 
the line. Tiny, barely-visible bubbles arc seen in most high flow rate runs, beginning 
some time after flow initiation, but these arc not large or numerous enough to 
compromise the experiment results. No erratic flow meter data were seen at any time. 

Samples of gas drawn from the ullage were analyzed with a mass spectrometer after the 
flight to detennine the constituents of the noncondensible gas fraction. Figure 23 shows 
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the results of the measurement, in percent by partial pressure. The measurement shows 
that the majority of the helium gas diffused out of the system during the ten months 
between the tank fill operation and the flight. During that time water vapor, nitrogen, 
and other components of air permeated into the tank, partially replacing the volume of 
lost helium. 



% of Noncondensibles 

Pressure. kPa 

Pressure, psi 

helium 

0.158 

0.65 

0.095 

water 

0.223 

0.92 

0.134 

nitrogen 

0.437 

1.80 

0.262 

oxygen 

0.155 

0.64 

0.093 

argon 

0.003 

0.01 

0.002 

carbon dioxide 

0.025 

0.11 

0.015 

Total: 

100 

4.14 

0.600 


Figure 23: Approximate Noncondensibles Concentrations During Flight 

Instrumentation. The lamps and cameras survived the flight environment and 
operated as planned, providing nearly four hours of real-time imagery. The quality of 
the photography is goxi and meets all the experiment requirements. Distortion and 
refraction limit visibility near the perimeter of the tank, though bubbles frequently 
provide a view of these regions by acting as a sphericd "fish-eye" lens. The timing 
LEDs are visible at the extreme edge of die picture. Both heaters are easily observable. 

The 3 megabytes of digital data were extracted and reduced in a multi-step process. 
Approximately 40 bytes were lost to transmission errors during this process. The data 
reduction software, which uses a mainframe con^uter to convert the binary data into 
formatted engineering-units data, has a bounds-checking feature that generates warnings 
whenever a value is above or below the range expected for a particular instrument. 

When this is encountered, the bad value is replaced by the previous value for that 
instrument and a warning is given. Approximately 500 warnings were generated from 
the STS-43 data, somewhat more than from ^und test data sets. These warnings, 
which are randomly distributed through the flight data and are not peculiar to a particular 
instrument, could be due to noise, transmission errors, or possibly cosmic radiation 
upset of tlffi memories (a known effect in the GAS radiation environment). In any case 
these represent a very small percentage of the data recorded (less than 6(X) bytes out of 
3,145,000 bytes). 

The pressure transducer signal was free of noise. Pressure oscillations observed in 
some phases of the experiment are believed to be real, consistent with heato* temperature 
swings, thruster firings, and flow-induced oscillations. As discussed previously, a 
2 kPa (0.3 psi) level shift was discovered in post-flight comparisons with three different 
calibrated gages; this shift is believed to have existed pritM" to the flight, so all flight data 
have been adjusted by subtracting this constant value from them. 

Evaluation of the flight data indicates that all seven thermistors in the fluid were 
operating no rmall y. The X-axis accelerometer's thermistor was consistent with the fluid 
thermistors to within 0.5* C at payload startup. The thermistor built into the computer 
reads consistently about 6-7*C higher than the accelerometer thermistor. This is 
acceptable since the computer temperature is fca: diagnostic information only. 

All three accelerometers performed within specifications, as nearly as can be determined. 
The rated accuracy is +/-0.12% in scale-factor and 4-/-0.008g bias (zero-offset), with 
one-year repeatibility in the bias of 4-/-0.(X)12g. The accelerometers were supplied with 
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factoiy-measured calibrations, which were confirmed by bench tests at Boeing prior to 
installation; however both of these measurements were made a year or mOTe prior to the 
flight The factory-supplied bias values for each unit (which is given as a function of 
temperature in the case of the X-axis unit) were subtracted from the data post-flight. 

The average acceleration values, calculated over each entire run after applying the 
factory-supplied conections, are plotted in figure 24. It may be noted that Aese time- 
averaged readings are not exactly zero in any of the three axes. This is almost entirely 
due to drifting of the bias values, or zero-offset, since the last calibration. All three units 
show a residual average reading of less than 0.(X)05g (absolute value). The amount of 
the drift is well within the one-year repeatibility specification for the instrument 
Therefore the absolute levels, or average values, of the accelerometer readings presented 
in this report do not indicate the true background g-level but represent residual 
instrument bias error. (However, as discussed in the next section, the change in these 
levels during the flight may be used to provide an indication of the background g level.) 

Note also that the average indicated values for the Y and Z axis units vary over the 
mission, as the payload warms up. The X axis unit exhibits the same trend before the 
temperature correction is applied. Therefore this trend is not real, but is a temperature 
effect since the Y and Z units are not thermally compensated. Note also a slight 
transient — ^within specifications — at the beginning of operation shortly after being turned 
on, due to instrument "warmup". 

Heater power and voltage measurements are consistent with ground tests and expected 
behavior. Heater A dissipates slightly more power than heater B (by approximately 0.5 
watts). The flow meter output is good and is consistent with ground tests and other 
data. Battery voltage at the end of the mission is consistent with predictions and 
previous behavior, thus indicating that the batteries performed as designed. 

The thermal environment was benign during flight operation. Payload temjjeratures 
recorded at the cortqjuter and at the accelerometers (near the tank at the lower end of the 
GAS caster) are plotted in figure 25. At tum-on 2 ^ indicated temperatures were 
approximately 19.5’C (67 °F), with the exception of the computer temperature which 
was 6*C higher. At the end of run 38, the accelerometers h^ warmed to less than 28°C 
(82'F), the computer indicated nearly 34*C (93*F), and the T7 trace indicates that the 
environment outside of the tank was less than 26*C (79‘F). The fluid mixed 
temperature in run 38 was 26.5'C (79.7*F), and the ma ximum fluid tenperature 
r^orded at any time was 34*C (93“F). These are significantly lower than the worst-case 
high temperature of 46°C (1 15*F) that was used to predict the ma ximum design pressure 
in the tank. 
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Figure 25: Payload Temperatures During Operation 


Data Summary 

Acceleration. The appendix contains plots of the complete set of digital data recorded 
during the STS-43 flight. Figures A-1 through A- 38 display the accelerometer ou^uts 
for runs 1 through 38, respectively. Because the absolute value of the data ranges over 
four orders of magnitude, each figure consists of two plots. The upper graph displays 
the data at a scale that includes the maximum range of the instruments. The lower graph 
plots the same data with an expanded scale so that smaller events and noise are visible. 
The same scales are used for each run so that they may be easily compared. The X-axis 
of all plots is the experiment time scale (the tinK elapsed since die beginning of run 1). 

Although the full-scale range of the accelerometers as read and recorded by the computer 
(uncorrected) is from -O.OlOg to -Kl.OlOg, the reduced data ranges from -0.014g to 
-t-0.(X)8g. This is because the raw output has been corrected in the post-flight data 
reduction by applying the factory calibrations for each instrument Again, note that none 
of the three channels' outputs are centered exactly on 0 g, due to the inherent bias drift 
of the accelerometers. 

The measured amplitude of the acceleration jitter is on the order of lO'^g peak-to-peak. 
Orbiter reaction control system (RCS) firings are visible in the data as nearly full-scale 
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spikes for the primary RCS thruster, and as 2 to 5 x IQ-^g spikes for the vernier 
thrusters. A firing of the orbit maneuvering system (OMS) thruster reaches full scale in 
both the X and Z axes during run 1 . 

Figure 26 shows the environment in run 3, a perhaps typical level of activity. The 
environment is quietest during the crew sleep period (approximately 240 minutes to 780 
minutes, TPCE-time). Figure 27 illustrates this, showing a significantly lower peak-to- 
peak noise level although vernier RCS firings are still seen. At other times of high 
activity, particularly during the TDRS/IUS deployment from 0 to 40 minutes, the 
environment is extremely noisy, as can be seen in figure 28. 
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Figure 26: Typical Acceleration Environment (Run 3) 
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Indicated Acceleration (G*s) 
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Figure 27: Sleep Period Acceleration Environment (Run 18) 
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Figure 28: Noi^ Acceleration Environment (Run 1) 


49 



Indicoted Accelerotion With Residual Bios (G's) 


Because the average acceleration level is very small in relation to the noise and to the 
range and accuracy of the accelerometers, it is not possible to directly determine the 
steady-state acceleration caused by drag and other background effects. It may howevCT 
be estimated in the X axis, or at least given an upper bound, by averaging the values for 
each run and observing the trend in relation to Goiter orientation, as plotted in figure 29 
(these values are Ae same as plotted in figure 24 for the X axis). If one assumes that the 
average acceleration due to atmospheric drag when the Qrbiter is oriented nose-first 
(+Xw attitude) is approximately equal and t^posite to that experienced when it flies tail- 
first (-Xw), then that value is one-half the difference between the average indicated 
levels for those two orientations. (Note however that time-averaging of the accelerations 
will tend to cancel out the random vibrations and crew motions, since they add no net 
inpulse to the vehicle, while a thruster firing which adds velocity to the vehicle will 
change the average value for the mn in which the firing occurs.) The approximate times 
during which the Orbiter was flown in the +Xw and -Xw attitudes, based tm the 
nominal mission timeline, is displayed on the graph for ccmparison. A horizontal line is 
drawn on the graph, somewhat arbitr^y, through runs 14 - 24, which occurred during 
the sleep period while in the -Xw attitude, and another is drawn through runs 25-31, 
which occurred after the Orbiter was maneuvered to the 4Xw attitude. The difference 
between these horizontal lines, 1.4 x 10-%, may be considered to be an upper bound of 
the step change in acceleration. Therefore the average acceleration due to drag is 
probably no more than one-half that value, or 7 x lfr%. 
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Tank Pressure (kPa) 


Pressure. Figures A-39 through A-76 in the appendix display the pressure, heat 
input, and mixer flow rate for runs 1 through 38, respectively. A typical example of the 
pressure history from runl 1 is shown in figure 30. 



Elopsed Time (Minutes) 

Figure 30: Typical Pressure History (Run 11, Mixer Flow 0.6 llmin) 

The pressure rise limit of 20.7 kPa (3 psi) is never reached during any heating stage, so 
each run's heating stage lasts the full ten minutes. Pressures are often seen to rise 
higher than they did in one-g testing, but the pressure rise is not snKX>th or monotonic as 
it was in one g. The pressure often exhibits little ot no change for several minutes, and 
then sharply rises to an unsteady plateau. In eight of the runs a spike is observed, 
where an increase of as much as 13.6 kPa (1.97 psi) occurs within two seconds. One 
of these spikes, recorded during run 13, is shown in figure 31. When they occur there 
is a simultaneous increase in the temperature recorded in the fluid nearest ^ active 
heater (T4 when heater A is active, and T5 when heater B is in use). In general, the 
longCT the period in which the pressure does not increase, the larger the subsequent 
spike. These spikes occurred when the video cameras were not operating, although 
smaller spikes are observed in several runs before the cameras are turned off. In run 5, 
for example, a small acceleration event approximately one minute into the heating stage 
causes the ullage to drift from the center of the tank to heater B, which is the active 
heater. A small pressure spike occurs at sixty seconds, approximately the same time that 
the ullage contacts the heater. 

When the mixer is turned on, ten minutes after the beginning of the run, the pressure 
drops rapidly as indicated in both figures 30 and 31. At low flow rates (such as figure 
30) the rate of pressure reduction is faster than observed with the same apparatus in 
noimal gravity tests. At all non-zero flow rates the pressure consistently collapses to 
near equlibrium, sometimes reaching a minimum value before beginning a slow increase 
as the pumps continue to run. In those runs having zero mixer flow, the pressure drops 
slowly (although faster than in one g), requiring five to fifteen minutes to approach 
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Tank Pressure (kPo) Pressure (kPo) 


equilibrium. An exception is run 1, in which the accelerations caused by the deployment 
of the TDRS^S payload (which is pushed out of the payload bay by springs) cause a 
sudden pressure r^uction. This self-mixing is shown in figure 32. The payload 
deployment produces a partial pressme reduction, and the OMS separation bum one 
minute later results in a more rapid pressure reduction. 



Figure 31 : Pressure Spike During Heating Stage (Run 13, Mixer Flow 2.8 l/rrUn) 



Figure 32: Self-Mixing Due To Accelerations (Run 1, Mixer Flow = 0) 
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Fluid Temperatures (C) 


Temperatures. Figures A-77 through A- 1 14 in the Appendix show the seven 
recorded fluid temperatures during runs 1 through 38, respectively. In each figure the 
upper graph plots the response of the four thermistors immersed in the fluid (Tl, T2, 
T4, and T5) plus that of Ae thermistor in the mixer loop (T7). The lower graph plots 
the two heater tenqjeratures (T3 and T6). The heater temperatures are plotted separately 
because their responses are slower than the fluid temperatures due to their thenr^ mass 
and therefore they do not provide an indication of the fluid mixing time. 

At the end of the heating stage the j>e^ fluid temperatures are typically 3-4’C (5-7°F) 
above equilibrium. When the mixer is turned on they converge within the 0.1 *C 
accuracy of the thermistors. The rate of convergence is a strong function of flow rate. 
Figure 33 is an example of the fluid temperature response to low-flow rate mixing (run 
1 1). The temperature response at a higher flow rate is shown in figure 34 (run 10). 

At the beginning of mixing in each ran, the temperature of the fluid in the mixing jet 
(T7) is approximately 1.0°C to 2.5*C (1.8°F to 4.5’F) below the equilibrium 
temperature of the tank contents. This is because when fluid is not circulating through 
the loop during the 25 minutes between mixing stages, it dissipates heat to the cooler 
environment Although the amount of fluid contained in the loop is small compared 
with the tank, this initially-cooler flow may have some effect on the pressure r^uction 
times. 



Figure 33: Fluid Temperature Response (Run 11, Mixer Flow = 0.6 Vmin) 
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Fluid Temperatures (C) 



Figure 34: Fluid Temperature Response (Run 10, Mixer Flow = 2.7 llmin) 

Video. The liquid/vapor initial conditions for each run were clearly capillary-force- 
dominated. The Bond number corresponding to the estimated 7x10*^ g-level is 0.01, 
for refrigerant 1 13 using the bubble radius of 13 cm (5 in). The ullage is generally a 
single spherical bubble that does not attach to the walls. One exception to this is that the 
accelerations during RCS firings occasionally force the bubble against a tank wall. 
Another exception is at the end of the heating phases, where the main ullage bubble is 
often attached to the heater, as will be discuss^ later. At other times the bubble drifts 
slowly, depending on the acceleration environment, but tqipearing to be repelled by 
protrusicHis such as the thermistor probes and the mixer nozzle. 

Figure 35 is a still photograph taken from the videott^je. It shows the ini^ conditions 
prior to heater tum-on at the beginning of run 9, typical of conditions during Ae crew 
sleep period for which the acceleration environment was most quiescent This was also 
during the -Xw ("tail-first") Obiter orientation. The image (and all subsequent video 
still photos) is oriented such that the mixer nozzle is to the right, heater A is at the left 
side, and heater B is at the top of the photo (although these components are not always 
visible in the photos). The uUage is visible toward the left of the photo ^ a smooth 
sphere. (Objects inside the tank are magnified by the lens effect of the liquid, and 
objects inside the cylindrical portion of the tank are also distorted and appear "wider" in 
this orientation; this is why the ullage does not appear as a perfect sphere.) 
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Figure 35: Initial Conditions: Quiescent Environment (Run 9) 

The -Xw attitude exists firom runs 4 through 25, approximately. In all but two of these 
runs the ullage bubble is initially located toward the nose of the Orbiter, as shown in 
figiOT 35, and is nearly motionless. In the remaining runs the initial orientation is 
variable. Thus it appears that the drag acceleration, ~7 x 10*7g, may be sufficient to give 
the fluid a prefened orientation. Although disturbances can cause the ullage to be 
displaced and to start drifting, it has a definite tendency to return to the end of the tank 
that is opposite the direction of flight 

Figure 36 shows the initial conditions during run 1, a time of greater acceleration noise. 
The effect of the disturbances on the surface can be seen as a significant amount of 
waviness and movement of the liquid/vapor interface. RCS thruster firings are 
sometimes seen to cause the ullage to begin drifting and, in a few cases, to flatten 
against a wall or heater plate. 
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Figure 36: Initial Conditions: Active Environment (Run 1) 

In most runs boiling is observed on the heaters during the heating phase. Boiling is not 
seen in runs 33 through 36, for which both heaters are operated with a total heat flux 
that is one-quarter that used in the other runs. Typically bubbles begin to be visible on 
the active heater approximately sixty seconds after the heater is turned on. The bubbles 
grow on both sides of the heater plate, remaining attached until they merge with other 
bubbles and detach or are otherwise dislodged. Figure 37 shows the conditions at the 
end of the heating phase in run 17, which used heater B. Several bubbles can be seen 
forming on the heater plate, and the main ullage bubble is also tangent to the heater. 
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Figure 37: Boiling ax End of Heating Stage (Run 17) 

Although only the first two minutes of each ten-minute heating stage are videotap^ the 
conditions at the end of the heating stage can be seen in the vidi^ recorded imme^ately 
prior to the initiation of mixing (as in figure 37). By this time, in the majority of runs 
exhibiting boiling, the main ullage bubble is attached or adjacent to the heater that is 
active for that run. The movement of the ullage fiom its initial location to the heater 
takes place during the non-videotaped portion of each run. 

In seven of the test runs a flow rate of zero is used to test the rate of pressure decay and 
re-equilibration in the absence of active mixing. The multiple bubbles present at the time 
the heater is turned off slowly re-combine into a single sphere by the time the following 
run is begvm (thirty minutes later). The coalescence of bubbles creates large visible 
ripples and occasional rapid movements in the ullage. This appears to produce a 
significant amount of self-mixing in the fluid, even during a quiescent acceleration 
environment. In the active acceleration environment of run 1, the interface distortion is 
more severe. 

The motion induced by the mixing jet is clearly visible in every run in which the jet is 
activated. Two main flow patterns are observed: ullage-penetrating and non- 
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penetrating. Both of these patterns are axisynmetric and are observed to form even 
when the initial location of the ullage is off-axis. 

At low flow rates the jet impinges on the bubble, creating a dimple or short geyser, but 
does not penetrate the surface. The bubble tends to be positioned and held at the 
opposite end of the tank, on its centerline, and the flow pattern is symmetrical and 
stable. Hgure 38 is representative of this nonpenetrating flow pattern. This photo is 
fiom run 1 1, which has an average flow rate of 0.59 1/min (0.16 gpm). At the lowest 
flow rate, 0.38 1/ min (0.10 gpm), the jet forms a barely-visible depression in the ullage. 
This is shown in figure 39, which is from run 25. 



Figure 38: Nonpenetrating Flow Pattern (Run 11, Mixer Flow = 0.6 l/min) 
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Figure 39: Lowest Flow Rate (Run 25, Mixer Flow = 0.4 Ihnin) 

At higher flow rates the jet penetrates the liquid/vapor interface, impinges on the 
opposite end of the tank, and recirculates along the tank walls. This results in a 
symmetrical annular bubble. An example of this is shown in figure 40. This photo is 
from run 12, which has an average flow rate of 1.54 1/min (0.41 gpm). This flow rate 
is at the low end of the ullage-penetrating range. At the high end the jet forms 
approximately the same pattern, but the hquid/vapor interface appears highly energetic 
compared with the glassy apjrearance in run 12. Figure 41 shows run 38, which has the 
highest flow rate reached, 3.35 1/min (0.89 gpm). During both of the runs with this 
flow rate the jet breaks the ullage up into a large number of small bubbles. 
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Figure 40: Penetrating Flow Pattern (Run 12, Mixer Flow — IS llndn) 
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Figure 41: Highest Flow Rate (Run 38, Mixer Flow = 33 llmin) 
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In every run having this symmetric penetrating flow, asymmetric flow is observed 
during some portion of the four minutes of video coverage. Often the flow begins as a 
symmetric pattern, degenerates into an asymmetric condition, and ftequently returns to 
the sy mme tric pattern. A run is classified herein as "penetrating" if at any time during 
the four minutes a symmetric, recirculating penetrating geyser is visible. One kind of 
asymmetric flow fir^uently seen at the higher flow rates is shown in figure 42 from run 
36, which has an average flow rate of 2.82 lAnin (0.75 gpm). Here the jet appears to 
remain on the axis but Ae ullage shape deviates from the near-toroidal shape to a two- 
lobed or saddle shape. The lobes move and in some cases stretch as far as the nozzle 
end of the tank. In another asymmetric flow pattern ftequently seen at moderate flow 
rates, the xiUage is pushed to one side of the tank but otherwise is shaped like the 
symmetric penetrating pattern. 



Figure 42: Asymmetric Flow Pattern (Run 36, Mixer Flow -2.8 Umin) 

At flow rates between those producing nonpenetrating and penetrating patterns, only 
asymmetric flow is observed. In this intermediate range the jet is genei^y deflected 
toward one side of the tank and the ullage is deflected toward the opposite side. 
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ANALYSIS AND DISCUSSION 


Flow Patterns 

The patterns observed in the videotape data may be con^aied with the four flow patterns 
i(fentified by Aydelott (refs. 13, 14, 15) that were discussed in the introduction. The 
low-velocity nonpenetrating regime is comparable to Pattern I, in which no geyser is 
formed. The bulge or dimple seen in the ullage at these flow rates does not qualify as a 
geyser — its width is greater than its height, and it does not significandy penetrate into 
the ullage. Flow Pattern n, in which the jet forms a long stable geyser diat extends into 
the ullage, is not seen, nor is Pattern HI, in which the geyser reaches the end of the tank 
and pools there without recirculating. At the startup of srane higher-flow rate runs, 
however, the jet penetrates the ullage at the initiation of mixing and forms a geyser, but 
the geyser quickly becomes either asymmetric (with the jet deflected to one side of the 
bubble) or reaches the end of the tank and forms the recirculating penetrating symmetric 
pattern. Thus it appears that Patterns II and m may not be steady-state flow patterns. 
Aydelott's drop tower tests were limited to 3.1 seconds duration, which may have 
prevented a stable flow pattern from developing fully. Moreover, the fluid initial 
conditions were symmetrical and undisturb^ which may have allowed an unstable 
symmetrical pattern to last the full three seconds of the test 

Flow Pattern IV is identical to the symmetric penetrating recirculating flow observed in 
TPCE. However, as discussed previously, in all runs that exhibit Pattern IV flow there 
are periods of asyinmetry. The flow frequently flip-flops between symmetric and 
asymmetric behavior, with no obvious acceleration events to cause dhe change. At these 
flow rates (>1.7 1/min) there does appear to be thorough mixing during both the 
symmetric and asymmetric periods— the jet sweeps all parts of the tank wall and is not 
blocked by the ullage from reaching any potential hot spots. 

As discussed in the introduction, the jet Weber number. We;, is a useful parameter for 
correlating low-g mixing performance. It is essentially the uirust of the jet divided by 
the resistance to penetration at the interface caused by surface tensirxi. Since the jet 
thrust apparently determines mixing effectiveness, and since Wej can be reduced to 
Aydelott’s flow characterization parameter and to several ullage-penetration criteria, this 
parameter can correlate both the flow patterns and the pressure control effectiveness. 

Use of Wej requires estimation of the jet diameter at the liquidtyapor interface, IX, 
accounting for spreading of the jet as it travels from the nozzle outlet to the interface. 

The Weber number calculations herein assume the following: 

1 . The jet spreads at 7° for the first 12.4 outlet diameters, and at 1 1* thereafter. 

2. The distance from the nozzle outlet to the interface, hb, is 12.7 cm (5 in). 

The first assumption is based on typical values of spread angle for submerged turbulent 
jets as given in reference 15. The second assumption is an approximation based on the 
following: for an 83% fill level, a spherical ullage has a diameter of 16.5 cm (6.5 in). 

If the ullage is located at the opposite end of the tank at the critical moment when the jet 
either penetrates it or is turned back, then the distance to the closest point of the sphere is 
12.7 cm (5.0 in). Using the nozzle outlet diameter of 1.0 cm (0.4 in) and the above 
spread angles, this results in a constant Dj of 4.94 cm (1.94 in). It is recognized that 
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this is a generalization and that D; may vary somewhat firom this approximation as the 
location and shape of the ullage changes, and that the spread angles are not necessarily 
correct when the jet impinges on an obstacle (the bubble). Rough scaling of the jet 
diameter from the video images does, however, confirm that this diameter is correct to 
within j^proximately 20% in both the symmetrical nonpenetrating and penetrating flow 
patterns. 

Figure 43 lists the resulting Weber numbers for each run having a non-zero flow rate, in 
ascending order, along widi the observed flow patterh. In this and subsequent figures, 
the nonpenetrating symmetrical pattern is identified by "Nonpen.", the peneuating 
symmetrical recirculating pattern is denoted by "Penetr.", and the asymmetric cases are 
identified by "Asym.". Nonpenetration is seen at Wej < 1.4, alAough at Wej = 1.4 the 
flow pattern is only marginally in this category. There is a gap in Weber number 
coverage between 1.4 and 3.1, which prevents an objective determination of the upper 
bound of the penetrating regime. However, since run 31 (Wej = 1.4) falls marginally in 
the nonpenetrating category, never appearing completely symmetric, while run 27 (Wej 
= 1.3) is symmetric and stable, it is believed that the upper bound of the symmetric non- 
penetrating flow pattern is at or near Wej = 1.4. 

Penetration (Aydelott's Pattern IV) is found for Wej ^ 4.8, altiiough there is an 
asymmetric ran at a higher Weber number. Run 28 (Wej = 5.9) is difficult to categorize 
but the flow pattern is never observed to be symmetric. At mixer startup, the ullage is 
attached to heater B where vigorous boiling is occuring. This asymmetric initial 
condition, caused by the geometry of the heat input, may have prevented prevented the 
jet firom forming the symmetric penetrating flow pattern during the four minutes of video 
coverage. Therefore, notwithstanding run 28, the lower boimd of the symmetric 
penetrating pattern is at or near Wej = 4.8. 

These Weber number boundaries may be compared with Aydelott's boundaries. If the 
low-g data from reference 15 are presented in terms of WebCT number rather th^ the 
flow characterization parameter F (and assuming that the jet Btmd number Boj is 
negligible), then Aydelott's reported flow patterns have the following ranges: Pattern I 
is firom Wej = 0.31 to 0.78, borderline Pattern VU. is seen ftom 0.39 to 1.11, Pattern II 
and n/En exists firom 1.16 to 1.59, Pattern HI is from 1.63 to 2.04, and Pattern IV is 
reported firom Wej = 1.33 up to 283. The upper bound of TPCE's nonpenetrating 
pattern (1.4) is only slightly higher than the upper bound of Aydelott's Pattern I/n 
(1.11). Both are somewhat higher, however, than Poth's criterion (ref. 2) which gives 
Wej^tical = 0-69. The bound^ between TPCE's asymmetric and penetrating flow 
patterns, Wej = 4.8, is higher than Aydelott's bounda^ between Patterns HI and IV 
(Wej = 1.33 to 2.04). 

The asymmetric range (1.4 < Wej < 4.8) produces irregular flow which, as will be seen, 
delivCTS inconsistent and sometimes less-effective mixing performance than either the 
penetrating or nonpenetrating regimes. At 3. 1 < Wej < 4.8, the jet is deflated to one 
side of the bubble and the return flow attaches to one side of the tank. This appears to 
be stable, since in this flow range the pattern is never seen to revert to a symmetric case. 
For 4.8 < Wej < 6.3, the flow alternates between the syrmnetric Pattern IV arid an 
asymmetric geometry, often taking the form of a saddle-shaped bubble. At higher 
Weber numbers (> 14.5) the ullage "saddle" becomes more pronounced, forming two 
peaks or arms, often with the flow also becoming symn^tric. Figure 44 is a sketch of 
the characteristics typical of some of the asymmetric flows encountered, as a function of 
Weber number. 
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Figure 43: Flow Pattern versus Flow Rate and Wej 



3.1<Wej<4.8 4.8 < Wej <6.3 14.5 < Wej < 22.6 

Figure 44 : Observed Behavior of Asymmetric Flows 
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Mixing Times 

When the mixer is turned, on the pressure drops at a rate that is dependent on the mixer 
flow. In some hi^-flow rate cases, the pressure reaches a minimum vtdue and then 
slowly increases. The Tninimiim pressure value reached at any time during the fifteen- 
minute mixing stage (after smoothing of the data to renoove high-fiequency fluctuations) 
is defined as Pniin- Figure 45 displays the pressure decay traces for four typical runs, 
normaUzed as (P - Pmin) / (Pinitiai - Pmin). This dimensionless pressure is plotted on 
both (45a) and logarithmic (45b) scales. The effect of flow rate on pressme 
collapse rate is evident in this comparison. The initial portion of the pressure ciure 
approximates an exponential decay in most of the runs, as indicated by the nearly 
straight lines on the logarithmic curves. 

The timpi required to reach the minimum pressure is difficult to extract fixon the dau 
because the flamess of the curve near the minimum and the amount of fluctuation in the 
data make the timft value highly susceptible to snu^ vpiatiais. However, this mixing 
timft is not of significant practical interesL In application, a pressure control system in a 
propellant tank that operates on demand would irobably not be designed to operate until 
tiie tniniminn pressure is reached. Because of diminishing returns as the pressure decay 
rate flattens out, it would probably be more efficient to operate a mixer until the excess 
pressure has been reduced by some fraction, or until it reaches the lower bound of a 
"dead-band" control pressure. 

Pressure collapse times and thermal equilibration times have often been defined as the 
rime required to bring the pressure or temperatures some percentage of the way to their 
equilibrium values. This procedure is used to derive pressure collapse times from 
TPCE’s data. (Some data smoothing was also used here to reduce the effects of random 
fluctuations and of resolution limits.) For instance, the 75% pressure collapse time, 
tp75%, is defined as the time elapsed from mixer tum-on until 


P(t) - Pmin =: 0 25 
Finitial " Pmin 


(13) 


This mixing time, along with the 50% collapse time and the 95% collapse time, was 
calculated 5>y a computer code) for all 38 runs. The 95% coU^se time proved to be an 
inconsistent measure of pressure control effectiveness due to numerical problenw caused 
by fluctuations in the curves. If one approximates the pressure response as having the 
form of an exponential decay, then the time constant of the pr^sirc decay can also be 
calculated as another measure of the pressure collapse rate which is potentially less 
sensitive to fluctuations in the data. This time constant, which is inversely proportional 
to the slope of the decay on a logarithmic plot (figure 45b), was also calculated for each 
run, by curve-fitting an exponential function to each decay. However this also proved 
to be prone to numerical problems and did not always represent the collapse rate well. 
The 75% pressure collapse time proved to be the least susceptible to numerical problems 
and is therefore presented herein as the most useful basis of comparison. 
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75Z Pressure Collapse Time (minutes) 


Figure 46 displays the 75% pressure collapse time, tp 7 s%, versus We; and flow pattern. 
In this and subsequent plots Wej is based on the average flow rate dumg the mixing 
stage, as CalciUat^ by a computer program. (Note however that in some rans the flow 
rate at the beginning of the mixing stage-where much of the pressure collapse occurs- 
vaiies by up to 5% from the average value over the entire stage.) Figure 47 plots the 
same information as figure 46 but on a logarithmic scale to better separate the high-flow 
rate mixing times. 



Figure 46: 75% Pressure Collapse Time Versus Weber Numer and Fbw Pattern 


68 



75% Pressure Collapse Time (minutes) 



Figure 47: 75% Pressure Collapse Time: Logarithmic Scale 

The higher flow rates — those that produce the penetrating symmetric flow pattern 
denoted on the plots by the square symbols — ^provide rapid pressure reduction, as 
expected. However, the low flow rates that do not result in ullage penetration (shown 
as circles) still cause the pressure to be reduced positively and consistently. Pressure 
collapse times are consistent for Weber numbers between 0.29 and 1.3. Although the 
nonpenetrating jet does not direcdy impinge on areas of the tank blocked by the ullage, it 
appears to impart significant motion to all surfaces of the bubble, which enhances heat 
and mass transfer well above that expected with conduction alone. The lone 
nonpenetrating data point that is well above the curve is run 31, which is at a borderline 
Weber number and displays asymmetric flow shortly after startup. Another data point 
that displays asymmetric flow (run 29) also falls si^ificantly above the rest of the 
points. A possible reason for this lack of repeatabihty in the asymmetric range is that 
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Jet Kinetic Energy Added to Reduce Pressure 75% (Joules) 


when the ullage deflects the jet to one side, the hot region may or may not be effectively 
cooled by the jet, depending on its location. Therefore asymmetric flow is potentially 
less consistent in its pressure control effectiveness. 

Pressme control efficiency — defined as the total jet kinetic energy expended in reducing 
the pressure 75% of the way to its equilibrium level — ^is plotted in figure 48. This is 
calculated by multiplying Ae jet power (or the kinetic energy added by the jet per unit 
time, (Ji/8)pVo^Do^ ) by the 75% pressure coUap^ time, Here an advantage of 
low-velocity mixing can be seen: the energy required at the low flow rates differs from 
that at the highest ones by approximately a factor of five. The longer time required to 
reduce the pressure at a low flow rate is more than conqwnsated for by the lower energy 
of the jet, because the kinetic energy addition rate scales with velocity to the third power. 
Again, runs 29 and 3 1 lie well above the trend of the other points. (Although the total 
amounts of energy added by TPCE's mixer are small, the mixer energy input could be 
significant for a large-scale cryogen tank operating for long durations in space.) 



Figure 48: Pressure Control Efficiency: Jet Kinetic Energy Added to 
Fluid to Reduce Pressure by 75% 
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Thermal equilibration times are calculated for each run in a manner similar to the 
calculation of the pressure collapse times. At any given time the maximum fluid 
temperature spread, AT, is defined as the highest measured temperature minus the 
lowest measured temperature (after again smoothing the data by time-averaging each 
data point with its adjacent points). TTiis calculation is limited to the four fluid 
temperatures in the tank: Tl, T2, T4, and T5. The thermal equilibration times are 
defined as the time required for this temperature spread to be reduced by a given 
percentage. The minimum value of AT measured in each run, or ATmin, is not zero 
because the O.rC accuracy of each thermistor results in an approximately 0.2*C spread 
between the temperatures when the fluid is well-mixed 

For example, the 75% thermal equilibration time, tx75%, is defined as the elapsed time at 
which 


AT(t)-ATmin (14) 

ATuiitiai - ATmin 

The thermal equilibration time constant was also calculated by fitting an exponential 
function to each run, as with the pressure, but again this method proved to be unreliable 
due to numerical problems. Therefore the 75% thermal equilibration time is presented as 
the best measure of bulk liquid mixing effectiveness. This is plotted as a fimction of 
Weber number and flow pattern in Figure 49. 

Since temperature is only measured at a few discrete locations, there is less repeatability 
in the thermal mixing times than in the pressure collapse times. Temperature, unlike 
pressure, is an intensive property and there is no assurance, or even likelihood that the 
four temperatures measu^ represent the maximum and minimum temperatures in the 
fluid. However the upper bound of the tr75% data envelope may be ccaisidered to 
represent the thermal equilibration effectiveness. This curve is similar in shape to the 
pressure curve. The pressure reduction times average approximately 30 percent less 
than the upper bound of the temperature reduction envelope, except at the highest flow 
rates where they are up to 90% less. 
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Figure 49: 75% Thermal Equilibration Time versus Weber Number and Flow Pattern 

Figure 50 Ksts the calculated mixing times for each of the 38 runs. These include the 
75% tfiamal equilibration time and the thermal equilibration time constant, plus the 
50%, 75%, and 95% pressiue collapse times and the pressure decay time constant. 

Also given is the time constant fitted to the temperature decay for the active heater in 
each run (i.e., if a particular run uses heater A, then the heater time constant given is 
based on fitting an exponential function to the T3 [heater A ten^CTatip] dam). The flow 
rates and Weber numbers given are average values for the entire mixing period. Figure 
5 1 lists the same data but sorted in ascending order of flow rate for ease in reviewing the 
trends and repeatability of the data. 
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Run 

Number 

Heater 


Weber 

Number 


! Thermal Mixing Times (Minutes): 

Pressure Reduction Times (Minutes): 

75% 

Time Const. 

Heater Const. 

50% 

75% 

95% 

Time Const. 

1 

B 

0.00 

0.00 

- 

4.18 

16.96 

4.30 

0.62* 

0.95* 

1.09* 

13.9* 

2 

B 

2.68 

14.51 

Penelr. 

0.18 

*•* 

0.62 

0.06 

0.09 

0.58 

0.09 

3 

B 

0.54 

0.59 

Nonpen. 

1.33 

9.23 

1.80 

0.54 

0.79 

2.32 

0.76 

4 

B 

1.53 

4.73 

Asym. 

0.02 

**• 

1.08 

0.17 

0.30 

0.47 

0.32 

5 

B 

2.72 

14.91 

Penetr. 

0.12 

0.35 

0.71 

0.04 

0.06 

0.09 

0.06 

6 

B 

0.00 

0.00 

- 

0.02 


2.74 

1.36 

2.89 

7,82 

2,03 

7 

B 

1.53 

4.74 

Asym. 

0.02 


1.13 

0.16 

0.29 

0.52 

0.21 

8 

B 

0.60 

0.73 

Nonpen. 

0.02 


1.61 

0.54 

0.72 

7,82 

1.20 

9 

A 

0.00 

0.00 

- 

7.28 

6.78 

3.23 

1.66 

3.44 

12.27 

2.66 

1 0 

A 

2.74 

15.16 

Penetr. 

0.57 

0.65 

0.56 

0.04 

0.08 

0.89 

0.07 

1 1 

A 

0.59 

0,71 

Nonpen. 

3.87 

4.86 

4.08 

0.33 

0.89 

6.22 

0.60 

1 2 

A 

1.54 

4.78 

Penetr. 

1.10 

1.22 

« « • 

0.09 

0.21 

3.57 

0.15 

13 

A 

2.78 

15.55 

Penetr. 

0.32 

0.56 

0.55 

0.04 

0.07 

0.46 

0.06 

1 4 

A 

0.00 

0.00 

- 

11.12 

22.57 

4.37 

1.10 

2.56 

6.71 

1.89 

15 

A 

1.53 

4.74 

Asym. 

0.82 

0.69 

0.74 

0.13 

0.37 

5.14 

0.27 

1 6 

A 

0.60 

0.72 

Nonpen. 

3.53 

3.64 

3.99 

0.43 

0.96 

7,19 

0.85 

17 

B 

2.78 

15.62 

Penetr. 

1.22 

0.27 

0.84 

0.07 

0.11 

0,30 

0.12 

1 8 

B 

0.00 

0.00 

- 

0.92 

0.80 

4.59 

1.06 

1.96 

4.56 

1.73 

19 

B 

1.58 

5.06 

Penelr. 

0.02 

*** 

0.80 

0.17 

0.30 

7,22 

0.25 

20 

B 

0.60 

0.73 

Nonpen. 

0.10 

*** 

1.59 

0.36 

1.05 

11.56 

0.71 

21 

A 

0.00 

0.00 

- 

7.08 

12.59 

5.46 

7.02 

6.13 

12.80 

42.7* 

22 

A 

2.84 

16.22 

Penelr. 

0.47 

0.57 

0.56 

0.04 

0.09 

3.47 

0.08 

23 

A 

0.62 

0.77 

Nonpen. 

2.22 

1.71 

5.77 

0.36 

1.18 

6.36 

0.81 

24 

A 

1.57 

4.96 

Penetr. 

0.55 

0.69 

0.81 

0.12 

0.29 

4.25 

0.22 

25 

B 

0.38 

0.29 

Nonpen. 

0.02 

•e* 

2.67 

0.79 

1.01 

9.47 

1.74 

26 

B 

1.24 

3.11 

Asym. 

0.15 

0.20 

1.39 

0.28 

0.47 

9.81 

0.54 

27 

B 

0.80 

1.29 

Nonpen. 

1.97 

16.25 

1.50 

0.29 

0.54 

8.58 

0.46 

28 

B 

1.71 

5.90 

Asym. 

0.02 

#** 

0.96 

0.09 

0.22 

0.53 

0.17 

29 

A 

1.24 

3.10 

Asym. 

0.32 

0.26 

1.15 

0.49 

1.54 

11.00 

1.23 

30 

A 

1.77 

6.30 

Penelr. 

0.68 

1.77 

0.83 

0.07 

0.12 

0.92 

0.09 

31 

A 

0.84 

1.44 

Nonpen. 

0.50 

0.49 

3.18 

0.50 

2.30 

11.44 

1.68 

32 

A 

0.38 

0.30 

Nonpen. 

0.77 

0.73 

3.38 

0.51 

0.87 

2.54 

0.83 

33 

AB 

0.64 

0.82 

Nonpen. 

1.43 

3.21 

1.47 

0.26* 

3.46* 

12.63* 

4.16* 

34 

AB 

1.54 

4.79 


0.07 

0.06 

1.19 

0.69* 

2.33* 

13.07* 

48.5* 

35 

AB 

0.00 

0.00 

- 

10.12 

39.80 

8.01 

0.09* 

0.09* 

0.79* 

0.29* 

36 

AB 

2.82 

16.08 

Penetr. 

0.03 


0.95 

0.18* 

0.22* 

0.77* 

0.35* 

37 

B 

3.34 

22.48 

Penetr. 

0.03 


0.65 

0.02 

0,04 

0.06 

0.04 

38 

A 

3.35 

22.64 

Penetr. 

0.27 

0.26 

0.41 

0.03 

0.06 

0.34 

0.06 


* values not meaningful (see text) *** value could not be fitted to data 


Figure 50: Summary of Mixing Times 

Note that certain pressure collapse values are flagged with asterisks in these tables. This 
indicates that the mixing time obtained by the data reduction program is not meaningful 
because of an insufficient initial pressure differentiaL In run 1 an acceleration 
disturbance causes complete self-mixing prior to the start of the mixing phase. In runs 
33 to 36, which use both heaters and produce one-quarter the heat flux, Ae initial 
pressure differential is small compared with the resolution and noise level in the data 
(see for example the pressure trace for run 34, figure A-72). These five runs are not 
included in the 75% pressure collapse time plots. In other cases, denoted by the 
code was unable to fit a time constant to the data. 
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Run 

Number 

Heater 




1 Thermal Mixing Times (Minutes): 

1 Pressure Reduction Times (Minutes): 

75% 

Time Const. 

Heater Const. 

50% 

75% 

95% 

Time Const. 

1 

B 

0.00 

0.00 

- 

4.18 

16.98 

4.30 

0.62* 

0.95* 

1.09* 

13.9* 

6 

B 

0.00 

0.00 

- 

0.02 

*** 

2.74 

1.36 

2.89 

7.82 

2.03 

9 

A 

0.00 

0.00 

- 

7.28 

6.78 

3.23 

1.66 

3.44 

12.27 

2.66 

14 

A 

0.00 

0.00 

- 

11.12 

22.57 

4.37 

1.10 

2.56 

6.71 

1.89 

1 8 

B 

0.00 

0.00 

- 

0.92 

0.80 

4.59 

1.06 

1.96 

4.56 

1.73 

21 

A 

0.00 

0.00 

- 

7.08 

12.59 

5.46 

7.02 

8.13 

12.80 

42.7* 

35 

AB 

0.00 

0.00 

- 

10.12 

39.80 

8.01 

0.09* 

0.09* 

0.79* 

0.29* 

25 

B 

0.38 

0.29 

Nonpen. 

0.02 


2.67 

0.79 

1.01 

9.47 

1.74 

32 

A 

0.38 

0.30 

Nonpen. 

0.77 

0.73 

3.38 

0.51 

0.87 

2.54 

0.83 

3 

B 

0.54 

0.59 

Nonpen. 

1.33 

9.23 

1.80 

0.54 

0.79 

2.32 

0.76 

1 1 

A 

0.59 

0.71 

Nonpen. 

3.87 

4.86 

4.08 

0.33 

0.89 

6.22 

0.60 

16 

A 

0.60 

0.72 

Nonpen. 

3.53 

3.64 

3.99 

0.43 

0.96 

7.19 

0.85 

8 

B 

0.60 

0.73 

Nonpen. 

0.02 

*** 

1.61 

0.54 

0.72 

7.82 

1.20 

20 

B 

0.60 

0.73 

Nonpen. 

0.10 

**• 

1.59 

0.36 

1.05 

11.56 

0.71 

23 

A 

0.62 

0.77 

Nonpen. 

2.22 

1.71 

5.77 

0.36 

1.18 

6.36 

0.81 

33 

AB 

0.64 

0.82 

Nonpen. 

1.43 

3.21 

1.47 

0.26* 

3.46* 

12.63* 

4.16* 

27 

B 

0.80 

1.29 

Nonpen. 

1.97 

16.25 

1.50 

0.29 

0.54 

8.58 

0.46 

31 

A 

0.84 

1,44 

Nonpen. 

0.50 

0.49 

3.18 

0.50 

2.30 

11,44 

1.68 

29 

A 

1.24 

3.10 

Asym. 

0.32 

0.26 

1.15 

0.49 

1.54 

11.00 

1.23 

26 

B 

1.24 

3.11 

Asym. 

0.15 

0.20 

1.39 

0.28 

0.47 

9.81 

0.54 

4 

B 

1.53 

4.73 

Asym. 

0.02 


1.08 

0.17 

0.30 

0.47 

0.32 

7 

B 

1.53 

4.74 

Asym. 

0.02 


1.13 

0.16 

0.29 

0.52 

0.21 

15 

A 

1,53 

4.74 

Asym. 

0.82 

0.69 

0.74 

0.13 

0.37 

5.14 

0.27 

12 

A 

1.54 

4,78 

Penetr. 

1.10 

1.22 

* • * 

0.09 

0.21 

3,57 

0,15 

34 

AB 

1.54 

4,79 

Asym. 

0.07 

0.06 

1.19 

0.69* 

2.33* 

13.07* 

48.5* 

24 

A 

1.57 

4.96 

Penetr. 

0.55 

0.89 

0.81 

0.12 

0.29 

4.25 

0.22 

1 9 

B 

1.58 

5.06 

Penetr. 

0.02 


0.80 

0.17 

0.30 

7.22 

0.25 

28 

B 

1.71 

5,90 

Asym. 

0.02 

**# 

0.96 

0.09 

0.22 

0.53 

0.17 

30 

A 

1.77 

6,30 

Penetr. 

0.68 

1.77 

0.83 

0.07 

0.12 

0.92 

0.09 

2 

B 

2.68 

14.51 

Penetr. 

0.18 


0.62 

0.06 

0.09 

0.58 

0.09 

5 

B 

2.72 

14.91 

Penetr. 

0,12 

0.35 

0.71 

0.04 

0.06 

0.09 

0.06 

10 

A 

2.74 

15,16 

Penetr. 

0.57 

0.65 

0.56 

0.04 

0.08 

0.89 

0,07 

13 

A 

2.70 

15.55 

Penetr. 

0.32 

0.56 

0.55 

0.04 

0.07 

0.46 

0.06 

17 

B 

2.78 

15.62 

Penetr. 

1.22 

0.27 

0.84 

0.07 

0.11 

0.30 

0.12 

36 

AB 

2.82 

16.08 

Penetr. 

0.03 


0.95 

0.18* 

0.22* 

0.77* 

0.35* 

22 

A 

2.84 

16.22 

Penetr. 

0.47 

0.57 

0.56 

0.04 

0.09 

3.47 

0.08 

37 

B 

3.34 

22,48 

Penetr. 

0.03 


0.65 

0.02 

0.04 

0.06 

0.04 

38 

A 

3.35 

22,64 

Penetr. 

0.27 

0.26 

0.41 

0.03 

0.08 

0.34 

0.06 


* values not meaningful (see text) *** value could not be fitted to data 


Figure 51: Summary of Mixing Times (Sorted by Flow Rate) 

Comparison With Mixing Models 

The dimensionless correlation of mixing times, tmixVoD(/D2tank> used by Wollen (ref. 

7) and others is applied to the 75% pressure coUapse data and plotted in figure 52. 
Wollen's experimental value of this parameter, a constant 3.46 at Reynolds numbers 
greater than 5(XX), is plotted in the same figure as a horizontal line. Although Wollen's 
correlation is based on equilibration within the liquid, rather than pressure i^uction, the 
pressure data agree approximately (to within a factor of three) except for two asymmetric 
flow cases (runs 29 and 31). Figure 53 shows Wollen's correlation versus the 75% 
thermal equilibration times. Here the dimensionless mixing parameter does not correlate 
the data well, and even increases the spread of the mixing times. 
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Dimensionless Temperature Collapse Time, Tt75*VoDo/Dtank**2 Dimensionless Pressure Collapse Time. Tp75*VoDo/Dtank**2 



Jet R^nolds Number (dimensionless) 


Figure 52: Pressure Collapse Times Compared with Wollen's Correlation 



Jet R^nolds Number (dimensionless) 


Figure 53: Thermal Equilibration Times Compared with Wollen's Correlation 
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Dimensionless Pressure Collapse Time, Tp75^Q/Vliquid 


Aydelott's (ref. 15) mixing time correlation (equations 5 and 6) is shown in figure 54 in 
conq)aiison with the 75% pressure collapse ^ta. It is compart with the 75% thermal 
equilibration time in figure 55. In both cases the data do not confirm the correlation, and 
the trends with increasing flow rate do not match the slope of the line. Note, however, 
that this correlation was based on the time required to distribute dye in the jet to the bulk 
liquid, and might not be expected to compare well with pressure collapse and thermal 
equilibration times. 



Figure 54: Pressure Collapse Times Compared with Aydelott's Correlation 
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Figure 55: Thermal Equilibration Times Compared with Aydelott's Correlation 

A simplified closed-form solution for the pressure collapse time was developed and first 
presented (in a slightly different form) in reference 20. This model makes the following 
assumptions: (1) the ullage gas expands isothermally, (2) the ullage consists only of 
propellant vapor, (3) the vapor is an ideal gas with constant specific heats and constant 
compressibility, and (4) there is no heat transfer from the ullage to the liquid. With a 
constant-temperature ullage. 


to 

dt m dt 


(15) 


where dnVdt is the mass change due to condensation. An energy balance on the liquid 
surface whose area is A yields an expression for the evaporation/condensation rate: 


to 

dt 


h A 
X 


(Tsat ■ Tiiquid) 


(16) 


where h is the heat transfer coefficient between the bulk liquid and the interface, X is the 
latent heat of evaporation, Tsat Is the saturation temperature corresponding to the 
instantaneous pressure, and Tiiq nirt is is temperature of the bulk liquid away from the 
interface. Using the Clapeyron equation, and assuming that the liquid density is much 
greater than the vapor density, one can ^proximate the dependence of Tsat on pressure 
as 


dTsat _ Tj - 

dP ^ 

^ Pi 


(17) 
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where Ti and pi are the initial temperature and density of the ullage gas. Since the final 
(or minimum) pressure is the saturation pressure corresponding to the bulk liquid 
temperature (which is assumed to be invariant), one can approximate that 

Tsat ■ ^liquid ~ ^ ^ (P " Pfinal) (18) 

Pi 

which when substituted into equations (16) and (15) yields 


dP__h A 
* X 


JiL. 
Pi ^ 


(P ■ Pfinal) 

pV 


(19) 


where P is the instantaneous pressure and V is the ullage volume. Since we are 
assuming that the ullage remains at a constant tempCTature during the pressure decay, 
(P/p) = (Pi/pi). The solution to equation (19), which is an exponential decay, is then 


P(0 ~ Pfinal 

Pi ■ Pfinal 


e-t/x 


where the time constant x is: 


( 20 ) 


Pi2 ^2 V 

“ h A Ti Pi 


( 21 ) 


The time required for (P - Pfinal)/(Pi - Pfind) to to 0.25 (which is the definition of 
the pressure collapse time used elsewhere in this report) is equal to 1.39x. Therefore 


1.39 pi2 X2 V 

tP75% - h A Ti Pi • 


( 22 ) 


Dominick (ref. 9) reported that the following expression from Jakob (ref. 21) adequately 
predicted the convection coefficient, h, between the bulk liquid and a flat liquid surface 
layer on which an axial submerged jet impinges: 


= (0.205) ReO-231 prjl/3 . 

Inserting equation (23) into equation (22) yields 


(23) 




(24) 


This equation suggests that mixing time is a strong function of the initial ullage mass 
times its heat of condensation, which determines tiie ener^ that must be dispersed into 
the liquid to accomplish the pressure reduction. The previous correlations do not 
account for these parameters. (If the ullage expansion is modelled as an isentropic 
process, rather than isothermal, then the mixing time equation also has the ratio of 
specific heats in the denominator. However, the isothermal assumption probably more 
closely approximates most situations, and is conservative.) 

Pressure collapse times predicted by equation (24) are calculated for TPCE using typical 
values for pi, Ti, and Pj. For A, the interface area acted on by the jet, the projected area 


78 


75% Pressure Coliopse Time (minutes) 


of a sphere having the same volume as the ullage is used, i.e., A = Jt ruUage^- This 
prediction is shown in figure 56 in compariscxi with the measured pressure collapse 
times. The closed-form prediction agrees well with the data, although it slightly 
undeipredicts pressure collapse times in the nonpenetrating flow regime. Given the 
generalizations and assumptions made in the mc^l, its agreement with the data points to 
tile potential value of a closed-form prediction such as this in estimating mixing times. 



Figure 56: Pressure Collapse Time Versus Closed-Form Prediction (Eq. 24) 
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Heater Plate Convection 

The heater plates are constructed with a significant amount of thermal mass to help meet 
STS flight s^ety requirements by limiting the maximum potential heater surface 
tfimpCTanires. Because of this thCTmal mass, the heater temperature decay does not 
occur at the same rate as the fluid teirqjerature equilibration. The heater temperature 
decay does, however, offer an opportunity to estimate the convective heat transfer under 
the influence of jet flow in a low-g two-phase fluid environment. This heat transfer data 
could be useful in future design studies for estimating the effectiveness of axial jets at 
cooling tank walls and other objects, as a function of the jet Reynolds number, flow 
pattern, and location of the area being cooled. 

A timft constant x is obtained for each run by fitting an e}qx>nential function to the^ 
temperature decay curve of the active heater. The ccmvectitHi coefficient h is obtained 
from this time constant by the relation x = hAAncp, where A is tiie total surface area of 
the heater plate (including both sides), m is its mass, and Cp is its specific heat. The 
coefficient is nondimensionalized by calculating the Nussdt number, hd^i, where d, the 
characteristic dimension, is the lengtii of the heater plate and ki is the liquid thermal 
conductivity. These Nusselt numbers are presented in figure 57 as a function of the jet 
Reynolds number at the outiet of the nozzle. The open symbols on the plot are heater A 
runs, and the filled symbols are heater B runs. 

The results indicate that heater A receives better cooling than heater B at high flow rates. 
This is not surprising because at these flow rates the jet penetrates the ullage and 
impinges directly on heater A, and heater B sees only a portion of the ret^ulating flow 
along the tank sidewall. However, at lower flow rates tlmt do not result in ullage 
penetration heater B experiences more cooUng from the jet than does heater A, which is 
blocked by the ullage bubble from direct jet impingement 
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Nusaelt Number Based on Heater Length (dimensionless) 
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Figure 57: Flow-Induced Cooling of Heater Plates 


Heating Phase Effects 

As discussed earlier, the main ullage bubble was observed to migrate to the active heater 
during nearly every heating period. This effect, also observed by other low-g 
experimenters, is due to Marangoni convection. As illustrated in figure 58, when a 
thermal gradient exists along a liquid/vapor interface there is a gradient in surface tension 
force that causes liquid to be pulled in the direction of lower temperature (which has 
higher surface tension). This results in motion of the bubble toward the heat source. 

The dimensionless Marangoni number is estimated to be on the order of 300 for TPCE, 
which is sufficient to drive this convection mechanism. The mechanism does not exist, 
thou gh, until the thermal gradient induced by the heater reaches the bubble. In most 
cases tMs probably takes several minutes, which is why the ullage migration is not 
observed in the two minutes of video at the beginning of each heating stage. 
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Figure 58: Marangoni Convection Mechanism 

The piessuie spikes seen in several of the heating phases, as discussed previously, are 
suspected of being related to this ullage migration. In at least two cases (runs S and 8) 
the ullage travels to the active heater early enough in the heating phase to be captured on 
the vidra because of an acceleration event. A small spike is seen in the pressure data at 
the same time that the ullage reaches the heater in these cases. It is not known if the 
larger spikes seen in other runs occur at the same time the ullage reaches the heater 
because there is no video coverage of those events. However this suggests a 
mechanism, as follows, to explain the spikes. 

When the ullage is not adjacent to the heat source in low-g, the heat that is added to the 
fluid cannot affect the pressure unless and until evapoiaticHi occurs. Boiling at the heater 
is suppressed in TPCE by the presence of noncondensible gases in the ullage and by the 
lack of nucleation sites. These factors mean that the liquid must be supeiheated above 
the refrigerant partial pressure existing in the ullage before evaporation can occur. Thus 
a significant amount of heat is added to the liquid surrounding the heater, creating an 
unstable situation. When the temperature gradient eventually reaches the ullage, 
Marangoni convection attracts the bubble toward the warmer fluid. As the liquid/v^por 
interface reaches the re^on of supeiheated liquid the excess energy results in sudden 
evaporation, or flash boiling. The pressure quickly falls back to a new equilibrium, 
however, as vapor condenses on the cold side of the bubble. 

This explanation is supported by the temperature and pressure data. The later in the 
heating phase that the pressure spike occurs, the greater the spike usually is, indicating 
that more energy had been stored in the liquid. An example is run 13, which uses heater 
A. Immediately before the pressure spike occurs the temperature of the liquid near 
heater A (T4) is 4.8’C above the bulk liquid temperature (die heater temperature is 14*C 
higher than the bulk liquid). A 4.8*C superheat conresponds to a 12.6 kPa (1.83 psi) 
vapor pressure differential. The pressure spike that occurs 7.5 minutes after heater turn- 
on has an amplitude of 13.6 kPa (1.97 psi). 

This phenomenon is not seen in normal gravity because buoyancy continuously draws 
the warmest liquid to the liquid/vapor interface, fri low gravity, however, the potential 
for pressure spikes may exist whenever localized heating occurs at a higher rate than the 
heat can be dissipated into the bulk liquid. Mixing would increase this heat dissipation 
beyond that expected with conduction alone and therefore should reduce the likelihood 
of these sudden pressure increases. 






CONCLUSIONS 


The Tank Pressure Control Experiment accomplished all of its objectives in flight on 
STS-43. The flow patterns produced by an axial jet mixer in a tank with an 83% fill 
fraction were characterized as a function of flow rate. The effectiveness of the mix ing 
was determined in terms of (1) the jet's ability to penetrate the ullage and reach all 
portions of the tank, (2) the time required to reduce the pressure in the tank, and (3) the 
time required to mix the fluid and equilibrate its teirperature. 

Two symmetric flow patterns were found: a nonpenetrating jet at jet Weber numbers 
less than approximately 1.4, and a penetrating, recirculating jet at Weber numbers of 4.8 
and above. Effective mixing was found with both flow patterns. At intermediate Weber 
numbers the flow was asymmetric and the mixing performance was sometimes less 
effective and was less repeatable. The highest flow rates, at Weber numbers near 23, 
resulted in breakup of the ullage into numerous small bubbles. 

Mixing times — both pressure reduction times and thermal equilibration times — were 
shortest at the highest Weber numbers, but the low end of the ullage-penetrating regime 
(Weber numbers between 4.8 and 6.3) still produced rapid, positive pressure collapse. 
This range could be used in applications requiring a rapid reaction of tank pressure. 
Based on the observed behavior at high flow rates, Weber numbers in excess of 
approximately 16 could result in undesired ullage breakup and fluid motions in large 
cryogenic propellant tanks. The nonj^netrating flow range (Weber numbers between 
0.3 and 1.4) produced a slower but still positive and repeatable pressure reduction. This 
range of Weber numbers might be used in applications requiring minimum energy 
dissipation. In this nonpenetrating range, the total kinetic energy dissipated to 
accomplish a 75% pressure reduction is approximately one-fif£ of that dissipated when 
the Weber number is in the ullage-penetrating range. 

Comparison of the results with prior work indicated that jet Weber number is a useful 
measure of both flow pattern and mixing effectiveness. This means that a target Weber 
number would be selected when designing a mixer to meet a given performance 
requirement One implication of this is that a mixCT can be designed to be more efficient 
by using a larger-diameter mixer nozzle which, at a given Weber number, would result 
in less lanetic energy addition. 

Previous mixing correlations did not accurately predict either the pressure reduction 
times or thermal equilibration times. However, a closed-form equation predicting the 
pressure reduction time, based on a simple thermodynamic model, approximately 
matched the results. This riKxlel has greater potenti^ for accurately predicting pressure 
collapse times than previous correlations because the effects of ullage volume and other 
important parameters are included. 

The results indicate that an active cryogenic tank pressure control system could be 
designed using a low- velocity axial jet mixer that would have an insignificant boiloff 
pen^ty associated with its kinetic energy addition. Low-energy mixing — either 
continuous or periodic — can provide positive and consistent p>ressure control. Mixing 
also offers the advantage of rnaking the fluid state tiK>re predictable and stable by 
reducing the thermal nonuniformities that can cause sud^n changes in pressure. 
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Figure A-95: Run 19 Temperatures 
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Figure A-96: Run 20 Temperatures 
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Figure A-97: Run 21 Temperatures 
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Figure A-98: Run 22 Temperatures 
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Figure A-99: Run 23 Temperatures 
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Figure A-106 
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Figure A-108: Run 32 Temperatures 
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Figure A-110: Run 34 Temperatures 
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Figure A-111: Run 35 Temperatures 
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Figure A-114: Run 38 Temperatures 
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Videotape S upplement 

A videocassette supplement to this document is available through the NASA Lewis 
Research Center videotape library. Much of the data returned by this experiment is best 
seen in this format as the still photographs reproduced in this report cannot convey the 
information that a motion picture can. The video supplement consists of two 2-hoi^ 
tapes containing the complete data set Each test run is identified in a preceeding utle 
slate by its run number, and a timer display is superimposed that corresponds with the 
time vines in the digital data. 

The videotape supplement is available on loan from NASA by sending the request 
shown below to: 

NASA Lewis Research Center 
Video Department (Mail Stop 5-2) 

21000 Brookpark Road 
Cleveland, OH 44135 


Date 

Please send one loan copy of (i^ Film Supplement to repKirt NASA CR-191012 

( ) Technical Film TF- ; (✓) Lewis Video 9M19_ 

in the following format (check one): 

( ) 3/4 in. U-Matic tape; (/) 1/2 in. VHS tape; ( ) 16 mm film 

Name of organization 

Street number 

City and State 2ip code 

Atm: Mr./Ms. — 

Title Telephone ( ) 




NASA Lewis Research Center 
Video Department (Mail Stop 5-2) 
21000 Brookpark Road 
Cleveland, OH 44135 
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